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An analytical model is developed to compare the effects of voids and debonds on 
the interfacial shear stresses between the adherends and the adhesive in simple lap 
joints. Since the adhesive material above the debond may undergo some extension 
{ either due to applied load or thermal expansion or both), a modified shear lag 
model , where the adhesive can take on extenslonal as well as shear deformation , is 
used in the analysis. The adherends take, on only axial loads and act as membranes . 
T wo coupled nondimensional differential equations are derived, and in general , five 
parameters govern the stress distribution in the overlap region. As expected, the 
major differences between the debond and the void occur for the stresses near the 
edge of the defect itself. Whether the defect is a debond or a void, is hardly discernible 
by the stresses at the overlap ends for central defect sizes up to the order of 70 
percent of the overlap region. If the defect occurs precisely at or very close to either 
end of the overlap , however , differences of the order of 20 percent in the peak 
stresses can be obtained. 


Introduction 

Adhesive bonded joints for composite or metal joining are 
being used in many structural applications. Among the most 
significant concerns regarding the structural design and reli- 
ability of the joints are the possible defects like debonds or 
voids which occur during manufacturing or service. These de- 
fects can severely reduce the bond strength. Their presence will 
increase the peak stress levels which occur at the joint ends 
and near the flaw itself. The joint may fail at the ends of the 
joint at the ultimate stress or it may fail under cyclic loading 
where local debonding near the flaw can grow. Also, the sub- 
sequent redistribution of stress, due to debonding, may lead 
to possible delamination in the composite adherend itself. The 
stress concentration, therefore, near a void or disbond is im- 
portant, and any thermal mismatch between the adherend and 
the adhesive will also contribute to the increased stress levels 
which occur. 

Past work, related to the present study, was performed by 
Erdogan and Ratwani (1971), Hart-Smith (1973, 1981), Kan 
and Ratwani (1983, 1986), and Rossettos and Zang (1993). 
The aforementioned work, dealing with voids or disbonds in 
the adhesive, uses a shear lag model, where the adherends take 
on only axial load and the adhesive takes only shear. This is 
appropriate in bonded joints which are designed so that the 
net load path does not produce bending. 

In the present work, which deals with simple lap joints, in 
order to properly compare the effects of debonds and voids, 
where the adhesive material above the debond may undergo 
some extension, it is necessary to allow the adhesive to take 
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on axial stress as well as shear. A modified shear lag model 
used by Rossettos and Shishesaz (1987) is adopted here for 
this purpose. As such, a quadratic distribution of axial dis- 
placement is assumed in the adhesive. Based on appropriate 
equilibrium, stress-strain and strain-displacement relations, the 
problem is reduced to two coupled second order differential 
equations for the axial loads in the adherends. A subsequent 
nondimensionalization of quantities in the equation leads to 
several parameters which are seen to govern the stress distri- 
bution in the joint. 

A structural mechanics rather than a continuum approach 
is used, where the loading mechanisms are restricted to net 
axiat and shear deformation in the components, and where a 
given quadratic displacement distribution is taken over the 
thickness of the adhesive layer. This avoids the corner sin- 
gularity at the overlap ends. Since the general solution of the 
structural model contains exponential terms, the steep stress 
gradients near the overlap ends and at the defect edges can be 
calculated accurately. It is also remarked that if transverse 
shear in the adherends is also included in a higher order analysis 
(Renton and Vinson, 1977) the peak shear stress values will 
occur very near the overlap ends, dropping sharply to zero at 
the ends themselves. The conclusions of the present paper, 
however, regarding these peak shear stresses, will change very 
little, if at all. 

Analysis 

The model consists of a simple lap joint as shown in Fig. 
l(tf, b) made of two plates, which take only axial loads, bonded 
by an adhesive layer. Plates 1 and 2 can be made of composite 
materials with orthotropic characteristics, although either ad- 
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herend could be specialized to an isotropic material with no 
change in the Formulation or the general nondimensional so- 
lutions. Let p, (jc) , Pi (x) and p 3 (x) be the resultant forces per 
unit width in adherend 1, adherend 2, and the adhesive re- 
spectively, while p 0 i$ the corresponding force applied to the 
adherends away from the joint. For any value of x we have 

Pi ( x ) + p 2 (x) +py(x) =p ( i (1) 

A quadratic displacement in the axial direction, W3 (a > , y), is 
assumed in the adhesive as 

ui(x, y) = u 0 (x) + [«2 W ~ u \ (*)J 

"0 

+ luz(x) +U\(x) -2z/ 0 (*)] (2) 

where u lt w 2 , and u 0 are the axial displacements in adherend 
1, adherend 2, and the center, ( y = 0), of the adhesive, re- 
spectively, and are functions of x. In Eq. (2), the origin of the 
local x-y system is at the center of the adhesive as shown in 
Fig. 1, so that u 3 at y = h 0 /2 and y = -h a /2 has the value 
w 2 and u l} respectively. In the adhesive, the lateral displacement 
in the x direction is taken uniform so that the shear stress is 
given by r(x, y) = G 0 du 2 /dy. On using Eq. (2), we obtain 

r(x, y) [w 2 U) “ u \ M]+^7T~ [“i (*) 
n 0 Hq 


+ u 2 (x)-2u 0 (x)] (3) 

By considering the equilibrium of each adherend, the fol- 
lowing equations are obtained 

^p-T 2 U)=0 -~T' +T i (•*) =0 (4a, b) 

dx dx 

where 

t 2 (x) = t ^x, J rj(x) = n y= -yj 

The shear stresses r 2 (x) and r,(x) are stresses in the adhesive 
just below adherend 2 and just above adherend 1 respectively, 


and can be calculated using Eq. (3). The net axial load in the 
adhesive is given by 


Pi(x) 


= ( T 
J h„ 


OxAx, y)dy 


(5) 


where <j , 3 is the axial stress in the adhesive and is a function 
of x and y. The axial stresses in adherends 1 and 2 are o x] = 
p x /h\ and o x2 = p 2 /7? 2 , respectively. Since the thicknesses h 0 , 
h x , hi are small compared to other structural dimensions, we 
assume the stresses o y \ = o y2 = o y2 = 0. We also assume 
generalized plane strain conditions to apply. These basic as- 
sumptions are also used by Kan and Ratwani (1983). The stress- 
strain relations are then given by 


b' X \“\ 

( 6a) 

*x2 “ r . ( 1 ” v zx2 v xz2 ) + a 2& T 

^x2^2 

(6b) 

y) i, . 

6^3 - (1 Vq) + Q-qA T 

(6c) 


where aj, a 2 , and <*o are coefficients of thermal expansion in 
adherend 1, adherend 2 and adhesive, respectively. The strain 
displacement relations are 




du i 
dx 1 




du 2 
dx' 


^3 = 


du 2 

dx 


(7) 


The principal equations for the problem can be derived as 
follows. If the derivative of Eqs. (4a, b) is taken, while ob- 
serving Eq. (3), derivatives of the displacements u 0 , «j, u 2 
appear. Equations (7) give these derivatives in terms of strains, 
and via Eqs. (6), (5), and (1), in terms of the loads p, and p 2 . 
In the process, it works out that the strain e x3 in the adhesive 
can be related to the derivatives of the displacements using 
Eqs. (6c) and (2), while the derivative of u Q can be eliminated 
in terms of Pi and p 2 using Eqs. (1), (5), (6c), and (7). In this 
way, the problem can be reduced to the solution of two coupled 
2nd order differential equations for the adherend loads p, and 
p 2 . The shear stresses at the interface between the adhesive and 
adherend are then given by the- derivatives of these loads as 
indicated in Eqs. (4a, b). 

Nondimensional quantities are next defined as follows: 

(Pi, P->) 

P( i 




(U u U 2 , U 0 ) = (Ml. MI. Mo) 

Pob 0 

_ dP L dp L ( h 0 \ 

=— -f = — T(x,y =±— ) (8) 

d f Po dx p 0 \ 2) 

where the quantities P, S represent nondimensional load, and 
shear stress, respectively, and will appear in what follows. With 
the definitions in Eqs. (8), the two differential equations are 
written in terms of the nondimensional loads Pi and P 2 as, 

^jr- Q(i + s,)P, - n p 2 = - fid - *,) (9a) 


d*P: 
d f 2 


■0(1 +5 2 )P 2 = -0(1-*) (9b) 


where 


0 = 3(1 


■"'(i 


. _2 Epho 1 
' 3 E tX hi 1 Mo 


5 2 = 


2 Epho 1 - u^iVxzz 

3 E x ih 2 1 — vl 


and 
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• i 4. ■.•■feS-f! 


"1 


2(2oq + g 2 “ 3<*o) Co// 0 A T 
3(1 -v 0 ) Po 

2(g 1 + 2a; — 3ap) Gq/iqAT 

3(1 -v 0 ) po 


The parameters 5, (/ = 1, 2) and Q represent material and 
geometric properties of the adherends and the adhesive. The 
parameters <£,- (/ = 1, 2), are thermal parameters which also 
include adhesive properties and the load p 0 . 

A solution to Eqs. (9) can be developed by means of an 
eigenvector expansion. As such the homogeneous form of Eqs. 
(9) can be written in the matrix form 

ir-LU = 0 (10) 

where U r = I P u P 2 I and L is a square matrix involving Q, 
5j, and <5 2 . A solution to Eqs. (10) is then assumed in the form 

U = He x{ (II) 

where U and H are vectors of second order. Substituting Eq. 
(II) into Eq. (10) results in the eigenvalue problem (L - \ 2 I)H 
= 0 which leads to eigenvalues, X?, with corresponding eigen- 
vectors, H\ For this problem, all the eigenvalues are real, and 
contain the parameters 5,, 8 2> and Q. The homogeneous so- 
lution can then be written as 

U = H m (Ci<? x ' 5 + C 2 e- X ' { ) + H ,!i ( Cj« X;5 + C 4 e ~ x -‘) (12) 

The particular solution, which is a constant, is easily obtained 
from Eqs. (9) and contains the parameters 5, and <J> h 
In the solution procedure, the adhesive is divided into three 
regions. Equations (9) apply 'to the first and third regions, 
which border the ends of the overlap. In the second or defect 
region, which falls between the other two, the governing equa- 
tions will differ from Eqs. (9). General solutions in all regions 
will give rise to constants of integration which are determined 
by appropriate boundary and continuity conditions. 

The general solution to Eqs. (9) is 

P|(?) = C|<? x ' f + C 2 e - x ' f + C,e x - { + C 4 e'^ ( + A t (13a) 
P 2 (£) =M' , (C,e x > f + C 2 <T x ' f ) 

+ M 2, (C 3 e X2f + C 4 e' X2f ) + Ai (136) 

where A , and A z are particular solutions. The nondimensiona! 
shear stresses in the interface between the adhesive and ad- 
herends as determined from Eqs. (4) and (8) are given by 


S i(«) = - 


dP i. 
*£’ 


S 2 (£) = 


dPi 

di 


(14) 


Void Model. The general solution, namely Eqs. (13), ap- 
plies to regions I and III, and gives rise to 8 independent 
constants of integration. The boundary and continuity con- 
ditions are now given as follows. In the notation for the loads, 
Pij, the first subscript / = 1,2 represents the adherend, while 
the second subscript j = I, II, III represents the region number. 

In region I (0 s £ ^ £ |) we have 
i 3 u(0)=I;P 21 (0) = 0;P II (fr) 

=Pm(fr);p 2 i(^D=p2n(fr) (is) 

In region II (void) (£i < £ < £ 2 ), since the shear stress on the 
adherends is zero, Eqs. (4a, 6 ) imply that the loads are constant 
so P 1M = C 5 and P 2 n = C 6 . In this region the load in the 
adhesive is zero so that Eq. (1) leads to 

Pm + Pm = Cj + C 6 - 1 (16) 

In region III (£ 2 ^ ^ I) 

Pm(? 2 “)-Pm.(? 2 + ); P 211 ( s 2 ~ ) = P 2111 (£ 2 * ); 

Pimd) = 0; P 2 „,(I)=I (17) 

There are ten unknown constants of integration but only 
nine conditions in Eqs. (15), (16), and (17). An additional 


condition can be obtained as follows. The nondimensional 
shear quantities 5, and S 2 can be written in terms of nondi- 
mensional displacements using Eqs. (3) and (8). When they 
are evaluated at £ 2 + an d £f they can be written as 

5 l ($ = ^ + ) = [4(/ 0 (£ 2 + )-3C/ J (jn-(/:(^)] (18a) 

52(£ = £2 + ) = -[4C/ 0 (f 2 + )-C/,(^)-3t/ 2 ({ 2 + )] (186) 

5 l (£=£r)-[4t/o(£r)-3c/ l (tf)-(/ 2 (fr)] cm 
s 2 (£=£n = -[4c/ 0 ({r)-t/i(£r)-3CA(fr)] <i8<o 

The changes in the inter facial shear stresses, as we go from 
£f to £ 2 + over the defect region, are given in terms of dis- 
placements as 

AS|=S,(£ = f 2 +)-Si(£ = £r) 

=4[(/ 0 (i 2 + )-c/ 0 (?r)]-3[c/ l (f 2 + )-t/ I (fr)] 

-[(/:(£: + )-(/:(£D] (19a) 

and 


A5 2 = S : ({ = f 2 + )-S 2 ($ = $r) 

= “4((/ 0 (^)-f/ 0 (?r)] + [^(£2 + )-c/,(fr)] 

+ 3 (t/ 2 (? 2 + )-t/:(ff)] 096) 
Also, the change in the displacements of the adherends, as we 
go from to £ 2 , can be derived by noting, for the void case, 
that e x[ and e x2 are constants; this gives 


U\ (£:) - U\ (£1) “ 


GgL 2 

Po^o 


(a, (.v 2 ) -!/,(*,)] 


G 0 L - 
Poho 


r x 2 


dx 


1 


E rl h t 


PoP\u+ot\i\T 


CoL 2 


Po^o 


(?:-?.) (20a) 


Uz(h)-U. i(f,) = 


GJL 2 

PA 


[w, (jt 2 ) — « t (JC, )] 


Gd? r 2 . 

= — — txidx 

Poho J., 


-[ 


I ~ VyaV.Tz2 

Exzhx 


P0P211 + <X 2 P.T 


CtyL- 
p 060 


($2-^.) (206) 


Using continuity of displacements in the adherends, adding 
Eqs. (19a) and (196) and using Eqs. (20) gives 


AS, + AS 2 = f — 2 1 

+ 2 2 poPin - 2a, A T+ 2a 2 AT 


Exihi 


GpL Z 

P o6 0 


(f 2 -?i) (2D 


On noting the general solution (Eqs. (13)), and the computation 
of shears via Eqs. (14), Eq. (21) can be written in terms of the 
ten constants of integration, and so provides the additional 
equation needed for their determination. 


Debond Model. An approach similar to the void model is 
used in the case of the debond, except now there is adhesive 
material above the debond in region II. Also, the quadratic 
displacement in this region is different from that in Eq. (2) 
which still applies to regions I and III. In region II, the quad- 
ratic displacement for u 3 must reflect the fact that the shear 
stress in the adhesive, at the debond (see Fig. 1(6)) location y 
- -h Q / 2, is zero (i.e., du 2 /dy = 0). Using this condition, the 
displacement in the adhesive in region II can be derived to give 

4 

wjU, y) = u 0 (x) +— (w 2 (.v) -u 0 (x)]y 

+ A [u 2 (x) -MoU)]>' J ( 22 ) 

3/i o 
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Wit6 this relation, the shear stress between the adhesive and 
adherend 2 is given in region II as 

'=('■'■ 2) ’it '" 1 '"" 1 <23) 

The governing equation for adherend 2 in region II is then 
given by 

^ ~ \ QP{ “ \ n (' + \ 6: ) Pz = _ \ n[ 1 - {U: _ ^ (24) 

Since the shear stress on adherend I is zero (due to debond), 
the load there is still constant so P U1 = C 5 . This is substituted 
into Eq. (24), and a solution for P 2!J is easily obtained, with 
two additional constants of integration. There is therefore a 

total of 11 unknown constants (Cj C 4 from region I, C 5 

. . . C 7 from region II and C% . . . Cn from region III) to be 
determined. The conditions in Eqs. (15) and (17) still apply, 
yielding eight equations. Two additional equations arise from 
the continuity of shear stress between adherend 2 and the 
adhesive at points £i and £ 2 , so 


s 2l (f={r)=s2i.(*=sr) 

(25a) 

S 2 mU £ = *2 - )=S 2I „(? = £ 2 + ) 

(25*) 


The eleventh required condition can be obtained in a manner 
similar to that of the void, except that in the debond model, 
in region II, 6 *i is constant but 6*2 is not. The change in the 
displacements of the adherends as we go from to £ 2 are now 
as follows. For adherend 1 , U i(£ 2 ) - £A(£ i) is the same as Eq. 
(20a) but for adherend 2 (using Eq. (6&) and the solution for 
Pin) get 

r j j r**2 

U 2 (i 2 )-U 2 (^)=~lu 2 (x 2 )-u 2 (x [ )]^— e x2 dx 

Po”0 Po*0 Jjr, 


1 - U- v2 v xz 2 

Ex!? >2 


j 

P 0 

\ 



6: 


(S:-s c ,)C 5 



+ 



C 6 


e 






\ 




+ a 2 Ar(J 2~ f i ) - 


I 


GpL- 

Po^o 


(26) 


Also, an equation, which is the counterpart of Eq. (21) for 
the void case, and which applies to the debond case can be 
derived to give 


ASj 4- AS 2 = 


-2p 0 


l-Vzx\Vxz\ 

Ex\h\ 


1 ~ Vzx2Vxz2 

E x2 h 2 


J 





With this additional condition, all eleven unknown constants 
of integration can be found, and the interfacial shear stress 
quantities can be calculated. The algebraic details can be found 
in (Lin, 1992). 


Results and Discussion 

As indicated by the nondimensional differential equations, 
the stress distribution in the overlap region is governed by the 
five parameters fl, 5 U 5 2 , <f> u 4> 2 , and the size and position of 
the defect. These parameters contain both geometric and ma- 
terial properties of the adherends and adhesive. Peak shear 
stresses are therefore dependent on the selection of geometry 
or materials or both. 

Typical distributions for the interfacial shear stresses S\ and 
S 2 are given in Fig. 2 for Q = 1000 and 4000. The defect is at 
the center of the overlap with size A£ = {2 - i\ = O- 4 ’ Values 
of the other parameters are as shown. It is clear that there is 
very little difference in the stresses at the overlap ends for the 
void and the debond. Any differences occur at the edge of the 
defect itself, with steep boundary layer behavior. In fact, these 
stresses can have differences of 100-150 percent. Note, for 



debond: void: no void: 



Fig. 2 Nondimensional shear stresses S,, S, versus { for two ft values, 
[t]: ft = 1000; [2]: ft = 4000; o, = o 2 = 0.01. 6, = b 7 = 0.01, = 0.3, 

$2 = 0.7. symmetric defect. 
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| Fig. 3 Nondimensional shear stresses S,, S* at special points versus 

1 A£ for three 0 values. [1]: fi = 1000, [2J: 0 = 4000, [3]: 0 = 8000; debond: 

I — , void: - ■ • = 0.01, = &2 = 0.01; symmetric defect. 


instance, that at £ = £, - 0.3, S\ decreases for the void case 
while increasing for the debond case, and the reverse is true 
for 5:. Both stresses show a tendency for steep edge effects. 
The structural mechanics (rather than continuum) approach 
in this paper does not address the type of singularity that exists 
at these edges. 

In Fig. 3, values for the shear stresses at the ends of the 
overlap and at the edges of the defect are plotted against center 
defect size, A£ . Again we see very little difference in the stresses 
at the overlap ends for the void and the debond, for a large 
range of defect sizes. Only when defect size is of the order of 
70 percent of the overlap length can we find a difference of 5 
percent in the case when ft = 1000. The stress values at £ [ and 
£ 2 (defect edges) are clearly substantially different for the de- 
bond and the void, for the full range of defect sizes. 

Corresponding results for the nonsymmetric defect are shown 
in Fig. 4. Here the defect size, A£ = 0.4, is fixed and the left 
defect edge, £j, moves from left to right (£1 = 0.1 to 0.55). 
The same general conclusions on the small differences in the 
stresses at the ends of the overlap for the void and debond 
hold. Notice however, that in Fig. 3 (symmetric case) the curves 
for a given defect, for different Q values, do not cross. But in 
the nonsymmetric case (Fig. 4) the curves for the stresses at 
the defect edges, do cross for different 0 values. 

A special case occurs when the defect is at either end of the 
overlap. If a void occurs at the overlap ends, the result is just 
a shorter overlap length and a larger value of physical shear 
stress. In the debond case this is not so' since the adhesive above 
the debond can take some axial load. The physical shear stress 
at £ = £ : is plotted against £ 2 , (the right edge of the defect) 
with £, = 0 in Fig. 5. While the peak stresses at £ = £ : in 
adherend 1 are larger for the debond than those for the cor- 




0 0.: t 0.4 0.6 0 0.2 0.4 <u 


Fig. 4 Nondimensional shear stresses S,, S* at special points versus 
for three 0 values. [1]: 0 = 1000, (21: ft = 4000, (3]: ft = 8000; debond: 
— , void: * • • *, o, = ©j = 0.01, 6, = S 2 = 0.01, nonsymmetric defect of 
size, a 0.4. 




$2 


Fig. 5 Shear stresses. r,/py>, r^py,, at special points versus $ 2 for =* 

0.0 = 4000, © t = © 2 = 0.01, = 5 2 » 0.01. debond: — , void: 
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fespondfhg void, they are smaller at the same value of £ : in 
adherend 2 for the debond. This is so, because the entire 
overlap length can share the shear stress. 
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On the Peak Shear Stresses in 
Adhesive Joints With Voids 


J. N. Rossetlos 3 ’ 5 and E. Zang 4,5 


Introduction 

Adhesive bonded joints have many advantages in terms of 
stress distribution, design flexibility, and simplicity of fabri- 
cation. It is known, however, that defects in the adhesive can 
severely reduce the bond strength. The presence of voids or 
disbond type flaws in the adhesive will increase the peak stress 
levels which occur at the joint ends and near the flaw itself. 

Past work related to the present study was done by Hart- 
Smith (1981) and Kan and Ratwani (1983). The work deals 
with voids or disbonds in the adhesive, is a structural mechanics 
rather than a continuum approach, and uses a shear lag model 
where the adherends take on only axial load and the adhesive 
takes only shear. This is appropriate in bonded joints which 
are designed so that the net load path does not produce bend- 
ing. 

In the present Note, it is shown how the stresses in a single 
lap joint with a void are completely characterized by two non- 
dimensional parameters. In particular, it is shown, how, for 
sufficiently large values of a parameter 8 (see Eq. 7), peak 
shear stresses can be essentially unaffected by relatively large 
central void sizes. Since the parameter 8 contains both geo- 
metric and material variables, it indicates explicitly how peak 
stresses can be controlled by selection of geometry and/or 
material. The influence of the location of the void is also 
indicated. It is remarked that if transverse shear is also included 
in a higher order analysis (Renton and Vinson, 1977), the peak 
shear stress values will occur very near the overlap ends-drop- 
ping sharply to zero at the ends (consistent with free adhesive 
ends). The conclusion of the present Note, however, will change 
very little if at all. 

Analysis 

The assumptions in the present model, with the geometry 
in Fig. 1, are similar to those by Kan and Ratwani (1983), so 
that the thickness variation of the stresses in the adherends is 
neglected. Also, in the z direction it is assumed that t iz - t^ 
= 0. In what follows, P\(x) and p 2 {x) are resultant forces 
per unit width in the x-direction in adherends 1 and 2, re- 
spectively, and p 0 is the resultant force in the adherends away 
from the joint region. The corresponding displacements are 
u x (x) and u 2 (x) . Since the adhesive is assumed to take on only 
shear, the shear stress there is 

r(x) = ( G/h 0 ){u 1 -u l ) t (1) 
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where G and h 0 are the shear modulus and thickness of the 
adhesive, respectively. Relations for the forces and stresses in 
the adherends are written as 

P\+Pi=Po\ <*\ x =P\/hu o-u^Pi/hi. (2) 

Because of the assumptions ( o iy = o 2y - 0), the stress-strain 
relations for isotropic adherends are given by 

ci* = (l -A){Po-Pi)/E\h\ 

*lr=d -»l)P 2 /E 2 h 2 (3) 

where E U E 2 are Young’s moduli and v u v 2 are Poisson’s ratios 
of the adherends. The strain-displacement relations are 

e \ x = du l /dx, tix - du 2 /dx (4) 

Equilibrium of an element of adherend 2 gives dp 2 /dx -t(x) 
= 0. The derivative of this equation gives 

d l p 2 /dx i - dr/dx = 0 (5) 

By then using Eqs. (l)-(4), Eq. (5) leads to an equation for 
p 2 . It can be written in nondimensional form as 


d 2 p 2 /d?-e 2 p 2 = -e 2 /(i +/?) 

(6) 

where 


e 2 = L 2 G ( 1 -»>?)( 1 + R ) /hoh 

(7) 

tf=(l -?!)£■,&,/( 1-,?)£ 2 /!j 

(8) 

and where the nondimensionalization is given by 


(F,. Pi) = (p 2 , P\)/p 0 \ i =x/L 


S = Lr/po ; (U 2 , U l ) = GL( u 2 , u , ) /p 0 h 0 . 

(9) 

By observing Eqs. (7) and (8), the parameter, 0, 

involves 


overlap length and adhesive/adherend geometric and material 
properties, while R measures the degree to which the adherends 
are different. Normalized void size is given by £ 2 - £1 = 
q 2 /L - a\/L. Note that if plate 2 is an orthotropic composite, 
then without change in the basic Eq. (6), R is given by 

R=(\-vi*viz)E x h l /(\- v })Et x h 2 . (10) 

In the solution procedure, the adhesive is divided into three 
regions. Equation (6) applies to the first and the third regions. 
In the void regions, r = 0, so Eq. (6) is replaced by d 2 P 2 /dt 2 
= 0. The solutions to these equations in the three regions lead 
to six constants of integration which are determined by bound- 
ary conditions and appropriate continuity conditions between 
regions. 

In region 1, (0 < £ < £,) 

J>2i(0) = 0; P 2 1 ( £ j ) = P 22 ( £ i ) * (Ha, b) 

In region 2, (void) (£] s £ ^ £ 2 ) 

Pn = constant (12) 

P23(t2)=P2,({,)+(t2-t,)F(P 2 2- 1/(1 +/?)), (13) 

where ( )' - d( )/d £. Equation (13) represents the change in 
the shear stress from £ j to £ 2 because of unequal displacements 
of the adherends over the void region (Zang, 1990). 

In region 3 (£ 2 :£ £ < 1), 

^22(f2)=/ > 23({2); />23(1)=1. (14<7 , b) 

Note that in the notation, P 2h the subscript, /, denotes the 
solution for P 2 in region /. 
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Fig. 2 Shear stress factor (SSF = ratio of peak shear stress to average 
shear stress) versus normalized void size, symmetric void; R 

= 1 



Fig. 3 Nondimensional shear stress, S(£), versus l for both a void lo- 
cated near the left end of the overlap and no void (continuous curves); 
R = 1 


Numerical Results and Discussion 

As an example, for identical adherends ( R = 1), the peak 
shear stress is shown to depend on both the location and size 
of £he void and the value of 0. For the larger values of 8 , such 
as $ = 10 or greater, a void located at the central portion of 
the overlap region with void sizes up to 70 percent of the overlap 
length (Fig. 2), will not affect the peak stress. For the same 
void located near one end of the overlap (Fig. 3), however, 
where a high stress gradient already exists, the peak stress can 
be increased by about 20 percentjfor 8 = 10). Note that even 
for the relatively large value of 8 = 10 the stress at the void 
edge, although not the peak stress, does increase markedly 
(i.e., by over 35 percent in Fig. 3). 

However, for smaller values of 8 , say 8 = 3, a void located 
at the central part of_the overlap will also affect the peak stress. 
For example, when 0 = 3, for a void one third the size of the 
overlap, an increase of 18 percent can be obtained for the peak 
stress, while the same void near one end (Fig. 3) can increase 
the stress by 39 percent. 

By observation of Eq. (6), it is seen that the variation of 6 
can be either due to changing the adhesive properties G and 
h 0t or the length of the overlap, Z,, or adherend properties. 
Yet, for a given value of 8 the non-dimensional shear stress 
distribution is the same. 


Similar trends (Zang, 1990) of the effect of the parameter 
0, also occur for values of R which differ from unity (i.e., 
different adherends such as composite/metal). In particular, 
a higher peak shear stress occurs at { = 1 , when adherend 2 
has a smaller stiffness than adherend 1, (i.e., R > 1). It turns 
out that as R increases from 1 to 10, the peak stress almost 
doubles in value. However, for larger R t the stress rapidly 
becomes asymptotic to a Fixed value, suggesting that adherend 
1 is essentially rigid relative to adherend 2. 
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First-Order Perturbation Analysis of 
Transient Interlaminar Thermal Stress 
in Composites 


Y. R. Wang 6,8,9 and T.-W. Chou 7,8 


The first-order perturbation analysis of the three-dimensional 
transient interlaminar thermal stresses of a symmetric com- 
posite laminate has been performed in this paper . Numerical 
results for a four-layer angle-ply laminate have shown that the 
first-order analysis was necessary for the solution of thick 
laminate (thickness-to-width ratio e S 0.02). 


Introduction 

This study is a supplement to the work of Wang and Chou 
(1989), in which the three-dimensional transient interlaminar 
thermal stresses of a symmetric composite laminate were ana- 
lyzed by a zeroth-order perturbation technique. Since the 
higher-order terms were neglected and it seems impossible to 
make an error analysis due to the mathematical difficulties, 
the present effort has been made to estimate the accuracy of 
the result obtained by Wang and Chou (1989) by including the 
first-order term in the solutions. 
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Nondestructive Evaluation of Adhesively Bonded Joints 
by Acousto-Ultrasonic Technique and Acoustic Emission 


Abstract 

Reliable applications of adhesively bonded joints 
require an effective nondestructive evaluation technique for 
their bond strength prediction. To properly evaluate factors 
affecting bond strength, effects of defects such as voids and 
disbonds on stress distribution in the overlap region must 
be understood. At the same time, in order to use acousto- 
ultrasonic (AU) technique to evaluate bond quality, the 
effect of these defects on dynamic response of single lap 
joints must be clear. The stress distribution in a single lap 
joint with and without defects (void or disbond) is analyzed. 
A 0 parameter which contains adherend and adhesive thick- 
ness and properties is introduced. It is shown for bonded 
joints with 0 £ 10, that a symmetric void or disbond in 
the middle of overlap up to the 70% of overlap length has 
negligible effect on bond strength. In contrast frequency 
response analyses by a finite element technique showed that 
the dynamic response is affected significantly by the pres- 
ence of voids or disbonds. These results have direct implica- 
tion in the interpretations of AU results. Through transmis- 
sion attenuation and a number of AU parameters for various 
specimens with and without defects are evaluated. It is 
found that although void and disbond have similar effects on 
bond strength (stress distribution), they have completely 
different effects on wave propagation characteristics. For 
steel-adhesive-steel specimens with voids, the attenuation 
changes are related to the bond strength. However, the 
attenuation changes for specimens with disbond are fairly 
constant over a disbond range. In order to incorporate the 
location of defects in AU parameters, a weighting function 
is introduced. Using an immersion system with focused 
transducers, a number of AU parameters are evaluated. It is 
found that by incorporating weighting functions in these 
parameters better sensitivities (AU parameters vs. bond 
strength) are achieved. 

Acoustic emission (AE) activities of steel-adhesive- 
steel specimens with 0 = 3.4 are monitored. Two different 
formats of energy vs. time have resulted, each correspond- 
ing to the perfect specimens or the specimens with void or 
disbond. The relative acoustic energy and the number of 
events at failure are found to be a means for predicting the 
bond strength. 
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H. Nayeb-Hashemi and J. N. Rossettos 

1. Introduction 

Over the past decade, much research effort has been 
expended and numerous test instruments have been devel- 
oped in seeking a solution to the problem of nondestruc- 
tive^ inspecting adhesively bonded joints (Williams and 
Zwicke, 1982; Rose, 1984; Dickstein et al., 1989, 1991, 
1992; Rose et al., 1983; Subramanian et aL, 1991; 
Williams et al. 9 1984). Several nondestructive evaluation 
(NDE) methods are recommended for the inspection of adhe- 
sively bonded joints, in addition to the well established 
methods for detecting localized flaws, voids, or delamma- 
tions. However, neither of these approaches nor more 
sophisticated ultrasonic methods (using frequency and time 
domain information), have been shown capable of abso- 
lutely assessing joint strength. This is due, in pan, to the 
fact that no single ultrasonic measurement is a unique func- 
tion of a single bond property; each ultrasonic measurement 
is sensitive to changes in several bond properties. However, 
the strength of the bond may or may not depend on these 
properties. A multidisciplinary approach is required that 
combines NDE, adhesive technology, and solid mechanics 
analyses to form a basis for a comprehensive quality assur- 
ance solution. 

Many attempts have been made to advance the state-of- 
the-art in flaw classification analysis by using techniques 
and concepts from pattern recognition. The underlying 
premise is built upon linear system analysis, which 
assumes that the ultrasonic input energy, as it varies with 
time, is modified by the bond structure. The theory then 
asserts, that if one has the system input and output (echo), 
then the modifying mechanism can be found and character- 
ized. Rose et al. (1983), using this concept, evaluated adhe- 
sively bonded aluminum-to-aluminum specimens by defin- 
ing an "a" parameter as the ratio of received signal to the 
transmitted signal in immersion scanning experiments. The 
feature "a" was defined and was of known value for bond 
quality discrimination (Meyer and Rose, 1974; Rose and 
Meyer, 1973). A low value of ”a H indicates good transmis- 
sion of the stress wave energy, while the high value of "a" 
indicates most of the energy is reflected at the interface. The 
values of "a" at several locations of the bond area were 
evaluated and the feature M (3 M based on the surface integration 
of the "a" parameter was defined as 
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P = -i T Jfa(x,y)dxdy 

4h " (1) 

= X 1 "-i w i a i (x ’ y)+R 

where h is the distance between scan points (square bond 
area), and Wj are the weighting factors. The surface integral 

values and log weighting function (Rose et al., 1983) from 
the center of each coupon were used to construct the feature 
data. These were related to the bond strength. The above 
equation and features can be modified by considering which 
area of adhesive is more responsible for the bond strength, 
and by proposing a weighting function which brings this 
factor into consideration. 

The initial signature techniques were based on the anal- 
ysis of ultrasonic signals in the time domain, and by mak- 
ing use of the arrival time and amplitudes. Gericke (1963) 
suggested that a source of additional discriminating charac- 
teristics might be found in the frequency spectrum of the 
returned echo. He indicated that if a short pulse, rich in 
spectral content, were used, the size, nature, and shape of 
flaws could be more readily determined from the frequency 
domain. He also proposed a method involving the use of 
two widely separated frequencies for inspection. The change 
of the returned pulse shape from one frequency range to the 
other could possibly yield information concerning the defect 
size, shape and orientation characteristics. Henneke and 
coworkers (Henneke et al., 1983; Duke et al., 1984, 1986; 
Sarrafzadeh-Khoee et al., 1986; Talreja et al., 1984; Govada 
et al., 1985) introduced several moments of the frequency 
spectrum as a means of damage evaluation in composite 
materials. A variety of additional signature techniques for 
ultrasonic examination both in time and frequency domains 
are introduced (Dickstein et al., 1990; Vary, 1987; 
Williams and Lee, 1987; Vary and Bowles, 1979; Vary and 
Lark, 1978; Nayeb-Hashemi et al., 1985). However, these 
parameters may not be effective when used in interrogation 
of the bond quality without understanding the effects of 
defects on both the received ultrasonic signature and bond 
strength. 

The higher order crossing (HOC) method for signal 
analysis, often called the zero-crossing or level-crossing 
method, has been recently developed and applied in bond 
quality assessments (Dickstein et al., 1989, 1992). Ultra- 
sonic echo signals were obtained from several specimens 
representing various adhesive or cohesive bond properties. 
HOC features were calculated from these signals and used to 
characterize the various conditions of the sample joints. 
However, no correlation between HOC parameters and bond 
strength were presented. 

Acoustic emission is another nondestructive evaluation 
technique which has been projected to have potential of 
predicting structural integrity. AE has developed rapidly 
over the last two decades as a nondestructive evaluation 


technique and as a tool for materials research. AE signals 
can take many forms depending upon the material and the 
failure mechanism in the material. AE signals from defects 
in composites and geological materials generally contain 
information at low frequencies, 0.5-100 kHz, since at- 
tenuation is relatively high due to the complexity of these 
materials. Signals of significance in metals and brittle 
materials contain information between 100 kHz and 2 MHz. 
In this range a good compromise is found for most testing 
applications because ambient noise is low. A number of 
techniques are employed to isolate valid signals form noise 
in the time domain. Times of arrival can be used to permit 
geometric elimination of obvious noise through gating, and 
acceptance of only those signals which arrive from a 
particular region of the structure. At the present time, one 
can find many presentations of AE data in the literature. In 
time domain, these include: AE (ringdown) counts, rms 
voltage, number of events, energy rate, rise time, event 
duration, amplitude distribution, and numerous others. 

It is also possible to analyze the frequency content of 
both burst and continuous types of AE. Relating such mea- 
surements to the source mechanism is an extremely com- 
plex problem, not only because of the specimen and trans- 
ducer resonances and frequency dependent absorption effects, 
but also because of effects caused by the methods of analy- 
sis employed. Indeed, some investigators would argue that 
the information required to distinguish between different AE 
sources, or to describe the nature of operation of particular 
sources, is simply not available in the frequency content of 
the AE signals, particularly using piezoelectric transducers. 
Other investigators believe that spectral analysis can be of 
great value and that with proper methods of averaging and 
smoothing of data, spectral analysis can be used to identify 
and discriminate AE sources (Heiple and Carpenter, 1983; 
Egle et al., 1981). 

As it was pointed out above, nondestructive evaluation 
of adhesively bonded joints requires comprehensive studies 
of the effects of defects, sizes, and their locations on the 
stress distribution and dynamic response of bonded joints. 
These analyses provide bases for more effective interpreta- 
tion of various AU data. The purpose of this paper is to 
present our efforts toward understanding factors affecting (1) 
the bond strength and dynamic response of adhesively 
bonded joints, and (2) bond strength prediction by conduct- 
ing various AU experiments. AE activities of steel-adhe- 
sive-steel, and graphite epoxy composite-adhesive-graphite 
epoxy composite specimens with various defects were also 
monitored during destructive tensile tests. The results of 
these investigations are presented in the following sections. 

2. Stress Distribution in Bonded Joints-Theory 

There are various bonded joint configurations. Most of 
the joint configurations are designed to transfer load in 
shear. Of all the various joints, the single lap joint is most 
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Fig. 1 Schematic diagrams of a step lap joint with a void, and a differential element 


commonly used in design. Also, theoretical and experimen- 
tal investigations of the bond strength are less complicated 
with the single lap joint. For this reason, in this paper the 
effects of voids, disbond region, poor cohesive strength, and 
poor adhesive strength on the strength of the single lap 
joints are studied both experimentally and theoretically. 

The strength of a given type of joint depends on the 
shear and peel stress distributions in the bonded area. These 
stresses depend on the adhesive and adherends' properties and 
geometries. Assuming the adhesive layer to be homogenous 
and free of defects Volkersen (1938) derived the shear stress 
distribution in single lap joints, using shear lag model as- 
sumptions. Goland and Reissner (1944) modified Volker- 
sen’s analysis by considering the bending moment and peel 
stresses in the theoretical analysis. Other researchers have 
also obtained the shear and peel stresses in the bond joints 
by Finite element analyses (Ishai et al., 1977; Adams and 
Peppiatt, 1974; Hart-Smith, 1985, 1987; Cooper and 
Sawyer, 1979; Allman, 1977). There have been very 
limited theoretical investigations on the effects of defects on 
the stress distribution and bond strength in single lap 
joints. Since the bond strength may or may not depend on 
the defects, the stress distribution in a single lap joint with 
a void is derived in this section, using shear lag model 
assumptions. The details of the analyses can be found in 
Rossettos and Zang (1993) and Hashemi and Rossettos 
(1990). Here, some aspects of these analyses are briefly 
described and important conclusions are elaborated. 

The two dimensional model consists of a simple lap 
joint as shown in Fig. 1, made of two plates bonded by a 


layer of adhesive. The problem is formulated under the fol- 
lowing assumptions: 

a) The thicknesses h 0 , h lf h 2 , H are small compared to the 

other dimensions of the structure so that the individual 
layers may be under generalized plane stress (i.e. o ly * <J 2 y 

= 0). 

b) The thickness variation of the stresses in the plates will 
be neglected under the usual assumption that the surface 
shear stress transmitted through the adhesive layer acts as a 
body force. 

c) In the z direction (i.e. transverse direction), plane strain 
conditions will prevail such that t z = t\ z - E2z = 0* 

Referring to Fig. 1, let p*(x) and p^x) be the resultant 

forces per unit width in the x-direction in plates 1 and 2, 
respectively, and U](x) and u 2 (x) be the displacements in 

the x direction. The elastic properties and thermal expansion 
coefficient of isotropic material 1 are donated by Ej , vj and 
ttj, and those for isotropic material 2, are E 2 , v 2 and a 2 * 
The shear modulus of the adhesive is given by G. The joint 
is divided into three regions. The following variables are 
defined as: 

xj = beginning of each region (x x = 0, x 2 = a lf 
X3 = a 2 in region 1, 2, 3, respectively), 
region 1: 0 < x < a^ 
region 2: aj < x < a 2 (void) 
region 3: a 2 < x < L. 

The shear stress distribution in the lap joint can be 
obtained by writing the equilibrium equation of an incre- 
mental element in plate 2 (see Fig. 1) as 


3 



dP2 

dx 


-t(x) = 0 


( 2 ) 


where x(x) is the shear stress and pj is the normal force per 
unit width in the adherend #2. Assuming the displacements 
of the top and bottom of the adhesive layer are U 2 (x) and 
Uj(x), and using strain-displacement relations of 


E ix = 


du 


and e 2x = 


du 


dx ‘ A dx 

equation (2) can be written in the form of 

” 'T* ( e 2x “ e U ) = 0 • 
dx h„ 


(3) 


Because of the assumptions (a ly = Ojy = the stress- 
strain relations in the adherends are 


1 — v? 

- -rr 1 [p ° " p 2 + a i AT • 

Eihj 

1-v? 

e 2x =— -^P2(x)+ a 2 AT. 

E 2 h 2 


(4) 

(5) 


Substituting equations (4) and (5) into equation (3), and 
presenting it in a nondimensionalized form results in. 


where 


£^_e2p 2 =^i+vj? 
d£ 2 l + R 

(6) 

£ = x/L 

(7) 

(P 2 .P,) = — (P2.Pl) 

Po 

(8) 

ft 1-vj E l h i 
1-vf E 2 h 2 

(9) 

02=iL-|-(l-v?)(l + R) 

hohi Ej 

(10) 

T 2 G 

V= (a 2 ai)AT. 

Po h o 

(11) 


Equation (6) holds for region 1 and 3; in the void region x = 
0, so 


^=0 


( 12 ) 


The solution of equations (6) and (12) along with proper 
boundary conditions and continuity conditions between 
regions yields the normal and interfacial shear stress distri- 
bution in the adherends (Rossettos and Zang, 1993; Nayeb- 
Hashemi and Rossettos, 1990). 


3. Theoretical Results of Void Effects on Shear 
Stress Distribution in a Single Lap Joint 


The geometric configuration of the step lap jointjs 
given in Fig. 1 and the three parameters R, 6 and t|/ , 
which appear in the governing equation (6), and influence 
the stress distribution, are given by equations (9), (10), and 
(11). The close observation of the aforementioned equations 
reveals that a change in 0 involves a change in geometry 
and material properties of both the adherends and the adhe- 
sive. An increase in 0 may be due to any number of possi- 
bilities, such as an increase in the length of the joint and/or 
an increase in the shear modulus of the adhesive, a decrease 
in the thicknesses of the adherend 1 and/or thickness of the 
adhesive layer, a decrease of the Young modulus of plate 1 
and finally an increase in R. 

The parameter, R, itself represents only the geometric 
and material properties of the adherends. The decrease of R 
involves the decrease of Y oung's modulus and thickness of 
plate 1 or the increase of the properties mentioned but 
applied to plate 2. For similar adherends, it is equal to 
unity. 

The parameter, \j t, can be seen as a thermal parameter. 
It involves the temperature difference between the two 
plates and their respective thermal coefficients of expansion. 
The parameter, ij?, increases if the thermal mismatch 
increases, where the two plates have a large difference in 
coefficient of thermal expansion. It also increases when the 
joint length and shear modulus of the adhesive increases. 
Finally, ij?, is inversely proportional to the thickness of 
the adhesive layer and the axial load. 

The effect of 0 on shear stress S(£) = (L/p 0 ) t distri- 
bution is shown in Fig. 2. For an adhesively bonded joint 
with the same adherend materials (R=l) the results show, 
that the higher 0 value causes very uniform (almost zero) 
stress over the bonded area with the peak shear stresses con- 
fined to the small region near the edge. Figure 3 shows the 
effect of a symmetric _void on the shear stress distribution 
for different values of 0 . Here again, the peak stress distri- 
bution is confined to the edge of the bonded joint. 

To evaluate the degree of stress variation as it is 
affected by a void and 0 , a shear stress factor, SSF, has 
been defined as the ratio of peak stress over the average 
stress. The SSF has been computed for a set of increasing 
symmetric void sizes. Figure 4 shows that the SSF is con- 
stant over most of the void size range (up to 70% of the 
overlap length) for the higher values of 0 . So it is expected 
that the failure load would be independent of a void for this 
range of void size. Effects of thermal mismatch between 
adherends are reported in detail in Rossettos et al. (1991). 
The stress distribution is significantly changed for adherends 
with a large difference in thermal expansion coefficient The 
effect of a disbond on the shear stress distribution is also 
investigated. Here again disbonds up to 70% of the overlap 
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Normalized Overlap Length £ 

Fig. 2 Normalized shear stress S(£) versus normalized 
overlap length £, for adhesive joints with identical 
adherends (R=l) for several values of 0 . 



Normalized Overlap Length ( 

Fig. 3 Normalized shear stress S(£) versus normalized 
overlap length for adhesive joints with identical 
adherends (R=l) and a center void, for several values of 0 . 

length were found to have negligible effect on the peak 
shear stress for 0 > 10. 

5. Void Effects on Dynamic Response of Single 
Lap Joints - Theory 

In order to understand the effects of voids on dynamic 
response of lap joints, we have studied the linear frequency 
response of an adhesive joint with a void, with particular 
attention given to the overlap region. This will be useful in 
the effective interpretation of ultrasonic data for the bond 
quality evaluation. 

In the analysis, a finite element model is used to repre- 
sent a simple adhesive lap joint. A harmonic force excita- 
tion is applied at one end of the joint near the overlap 



Normalized Void Size 

Fig. 4 Shear stress factor ipeak^av versus normalized void 
size (£2 * £i)» in adhesive joints with identical adherends, 
and a center void, for several values of 0 . 

region, and a displacement response at several output points 
is calculated over a frequency range. The modal analysis and 
frequency response calculation is performed for a joint with 
and without a void in the adhesive to determine changes in 
frequency response patterns caused by the void. The results 
indicate clear differences in the response patterns over a fre- 
quency range which covers all natural frequencies of the 
finite element model. The results presented here focus on 
the higher frequencies. 

Beam elements are used for the adherends and lateral 
stiffness elements are used for the adhesive. Voids are cre- 
ated by removal of lateral stiffness elements. The model is 
simple and will indicate the important features to be 
expected. The NASTRAN code has been used. Material 
properties are for similar aluminum alloy adherends (E = 69 
GPa and adherend thickness of 3.17 mm) and Hysol 
EA9689 adhesive (E = 2.2 GPa and thickness of 0.13 mm). 
The overlap was 25.4 mm x 25.4 mm and symmetric and 
unsymmetric voids were introduced in the overlap and 
dynamic responses for identical input harmonic force were 
evaluated. 

The system equations formed for frequency response 
analysis include mass, damping and stiffness matrices and 
the system load vector. The equations may be written in 
matrix notations as follows: 

MU(f)+ CU(f)+ KU(f) = F(f) (13) 

where f is the discrete frequency value. 

The harmonic excitation is frequency dependent and the 
solution is obtained for the desired discrete frequencies. Since 
all nodal motion is assumed to be steady state, the ve- 
locities and accelerations are related to the displacements by 






FREQUENCY (MHz) 


Fig. 5 Finite element model and frequency responses at 
node 12 for aluminum-adhesive-aluminum specimens with 
no void, and 6.35 mm void. Void starts at 3.18 mm from 
left end of 25.4 mm overlap. 



Fig. 6 Finite element model and frequency responses at 
node 10 for aluminum-adhesive-aluminum specimens with 
5.08 mm and 10.16 mm void. Void starts at 2.54 mm from 
left end of overlap. 

U = -(2irf) 2 U, U = i(2ref)U (14) 

The frequency response equation then becomes 

[-(2itf) 2 M + i(2jif)C + K]U(f ) = F(f) (15) 

Equation (15) represents a system of equations and the solu- 
tion for U(f) can be found for every frequency point. 
Proportional damping is used so that 

C = cjK + C2M. (16) 

In the present study cj = 0 and C] is given very small val- 
ues (of the order of 0.001). Typical frequency response 
curves are given in Figs. 5 and 6. These are calculated for 
the two configurations indicated in the figures. For a given 
frequency, the ratio of the displacement with a void to the 



Fig. 7 Ratio of displacements at nodes 8 and 12 for 
specimens with and without void (A/A^ versus void size at 

frequency of 0.39 MHz. Void starts at 3.17 mm from left 
end of overlap, and void size is given as percent of overlap 
length of 25.4 mm. 



Fig. 8 Ratio of displacements at several locations for 
specimens with and without void (A/A 0 ) versus void size at 
frequency of 0.49 MHz. Void starts at 2.5 mm from left end 
of the overlap and size is given as percent of overlap length 
of 25.4 mm. 

displacement without a void is plotted against void size in 
Figs. 7 and 8. It should be borne in mind that for different 
void sizes, the structural configuration changes and the nat- 
ural or resonant frequencies will change, so that for a fixed 
frequency these curves need not be monotonic. 

The results of this analysis indicate that, although there 
was little change in the peak shear stress in adhesively 
bonded joints with voids up to 70% of the overlap length 
(0 £ 10), the dynamic response may be affected by the 
presence of voids significantly. These results have direct 
implication in the non-destructive evaluation of the adhe- 
sively bonded joints. 
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Fig. 9 System for attenuation and AU parameters measurements. 


6. Experimental Investigations 

Adhesively bonded joint specimens were prepared by 
using AISI 1018 cold rolled steel, unidirectional graphite 
epoxy composite material, and aluminum 6061-T6 as 
adherends. Identical adherends were joined together in a sin- 
gle lap joint configuration using Hysol EA9689 epoxy film 
of 0.13 mm thickness. For joints with steel as adherends, 
the adherend surfaces were prepared according to the ASTM 
standard D6251-79. For joints with aluminum as adherends, 
the adherend surfaces were either sanded and cleaned with 
acetone or were just cleaned with acetone prior to bonding. 
The composite adherends were slightly sanded in their over- 
lap areas, in order to remove the residual mold release prior 
to their bonding. For metal adherends the overlap area was 
25.4 mm x 25.4 mm and adherend thickness was 3.17 mm. 
Joints with graphite epoxy composite as adherend, had over- 
lap dimensions of 50.8 mm x 25.4 mm and the thickness 
of the composite adherend was 1.07 mm. Various defects 
such as voids, disbonds were introduced in the overlap area 
in order to change the bond strength. Voids of different sizes 
were introduced by cutting the adhesive film. The disbond 
was created by spraying mold release agent over the disbond 
area. 

Various ultrasonic setups and methodologies were used 
to measure longitudinal pulse-echo attenuation, through 
transmission longitudinal and shear waves attenuation, 
using both direct contact and immersion techniques, and 
acousto-ultrasonic parameters. In the attenuation experi- 
ments, narrow-band Panametrics transducers with center fre- 


quencies of 1 and 2 MHz were used both as transmitting and 
receiving transducers. Both narrow-band signal (pulse oscil- 
lator) and broad-band pulse were used in these experiments. 
In the AU experiments, narrow-band Panametric transducer 
with center frequency of 1 MHz and an FC-500 AET trans- 
ducer with a flat sensitivity in the frequency range of 100 
kHz-2 MHz, were used as transmitting and receiving trans- 
ducers, respectively. The experiments were conducted by 
using a broad-band pulse with the center frequency of 1 
MHz. The received signals were digitized in a Nicolet digi- 
tal oscilloscope and stored in a PC computer and later ana- 
lyzed in a SUN work station. Fig. 9. The signals were ana- 
lyzed in both time and frequency domains. 

In the time domain, attenuation, stress wave factor 
(SWF) and AU parameter (AUP) were measured for each 
specimen. Assuming the impedance between transmitting 
and receiving transducers and adherends are Fi(o>) and 
F 2 (oa), the amplitude or peak amplitude of the received sig- 
nal from a perfect and a defective specimens for an input 
pulse oscillator or pulse can be expresses as 

(Ar)p = F,(o))F 2 ( 0 ))A l exp(-apt) (17) 

(Ar) d = F,((o)F 2 (u>)A t exp(-a d t) (18) 

where A t is the transmitted signal amplitude, Op and a d are 
the attenuation in perfect and defective specimens, t is tht 
bond thickness, and (Ar)p and (Ar) d are the amplitude of the 
received signals from perfect and defective specimens at 
angular frequency to. The change in attenuation can be 
obtained from equations (17) and (18) as 
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Fig. 10 Schematic diagram of the AE data acquisition system. 


a J -a p = 4a=i(n[i^j (19) 

Stress wave propagation efficiency was evaluated by 
evaluating SWF values (Vary, 1987; Vary and Bowles, 
1979) at several threshold levels. The SWF is defined as the 
number of times a signal passes a threshold level. AUP 
were also measured for all specimens. This is defined as 

AUP = ^(A, -Th) (20) 

i=l 

where Aj is amplitude of the signal passing threshold level 
of Th, and n is equal to the SWF. Data were analyzed at 
different threshold levels including one which was set just 
above the noise level and similar conclusions were drawn. 

Upon completion of AU experiments, specimens were 
broken in an Instron testing machine by applying tensile 
load at the rate of 22 N/s. AE activities of stcel-adhesive- 
steel and graphite epoxy composite-adhesive-graphite epoxy 
composite were monitored. Here the results of AE activities 
of the steel-adhesive-steel specimens will be discussed. 
Since the failure load in these specimens were much less 
than the load required for yielding adherends, all received 
activities were related to the bond failure. However, for 
graphite epoxy composite-adhesive-graphite epoxy compos- 
ite specimens, bond failures were sometimes accompanied 


by fiber fracture and delamination. For this reason the corre- 
lation of these data with the bond strength was not clear. 

For steel-adhesive-steel specimens, two AE transducers, 
each with a center frequency of 175 kHz (AC 175L), were 
utilized to detect the events emitted from the specimen at a 
distance of 7.6 cm apart. Each sensor was coupled with the 
workpiece through B-type Panametrics couplant and held in 
position with the aid of 4 stiff springs as shown in Fig. 10. 
The data were gathered and analyzed using an AET 5500 
system. Data collection was focused on the lap joint por- 
tion only and all the other AE activities were thus discarded. 
For the lap joint, AE event locations were established using 
the difference in the arrival times of the signals received by 
the two sensors. Some of the parameters measured during 
the tensile tests were as follows; peak amplitude in deci- 
bels, event duration and rise time in (is, ringdown counts, 
slope, and AE energy. Here, the slope has been defined as 
peak amplitude/rise time while AE energy was defined as 
10*log(event duration) + peak amplitude (in dB). It is no- 
table that rms voltage has also been used as an indicator of 
the relative amount of AE energy by some researchers. 
Pencil-lead fracture method was used in order to calibrate the 
system prior to the actual run. Calibration eliminated the 
need to obtain the velocity of emitted waves. The detected 
signals were initially preamplified 60 dB to a total system 
gain of 80 dB. A floating threshold of 0.5 V was chosen to 
eliminate the background noise. 
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Fig. 1 1 Changes of longitudinal wave attenuation versus 
bond strength, for a pulse oscillator with the center 
frequency of 2 MHz, propagating through the lap thickness 
in steel-adhesive-steel specimens with symmetric void and 
disbond in the middle of overlap. 



Fig. 12 Changes of longitudinal wave attenuation versus 
bond strength for a pulse oscillator with the center fre- 
quency of 2 MHz, propagating through the lap thickness in 
steel-adhesive-steel specimens with symmetric void in the 
middle of overlap. 

7. Experimental Results and Discussion 

For the steel bonded joint specimens the value of 0 
from equation (10) is found to be 3.4. For this value of 0 , 
the peak shear stress is affected by introduction of void and 
disbond in the overlap. It was suspected that the through 
transmission attenuation would be also affected in a similar 
way by the presence of void or disbond. Using both 
Panametric longitudinal and shear transducers with nominal 
element diameter of 19 mm, the attenuation changes for the 
entire overlap length (25.4 mm x 25.4 mm) were evaluated 
at frequency of 1 and 2 MHz. Transducers were coupled to 
the specimens by using either low or high viscosity 
Panametric couplant (low viscosity couplant when using 



Time, sec 

Fig. 13 Typical load and ringdown counts versus time 
(interval size *= 2 s.) 

the longitudinal transducers and high viscosity couplant 
when using the shear transducers). A load of 90 N was 
applied to each transducer by using a spring fixture device. 
Figure 1 1 shows the changes in longitudinal wave attenua- 
tion vs. bond strength. Similar pattern was also observed 
when plotting shear wave attenuation vs. bond strength. At 
first glance the data suggest no correlation between bond 
strength and attenuation change. However, closer examina- 
tion of the data shows that, although disbonds and voids 
have similar effect on the bond strength, they have com- 
pletely different effects on wave propagation characteristics. 
For specimens with a void, the transmitted wave is reflected 
significantly at the adherend void interface. A larger void 
results in more reflection and thus higher attenuation. 
However, the disbond is a weak bond between adherend and 
adhesive. The transmitted wave is not significantly affected 
by the presence of this region in the overlap. Furthermore, 
the attenuation change may not be very sensitive to the dis- 
bond length. Separating the data for specimens with voids 
and specimens with disbonds, and plotting attenuation 
change versus bond strength results in an excellent correla- 
tion between attenuation change and bond strength for spec- 
imens with voids; see Fig. 12. There was no correlation 
between attenuation change and bond strength for specimens 
with disbond, and the attenuation change was fairly constant 
over the disbond of 3 to 19 mm range. The correlations 
between the SWF and AUP, and the bond strength for these 
specimens were not clear. This may be explained by consid- 
ering the dynamic response of the systems for an input 
wave at one edge and received wave at the other edge, as 
shown in Figs. 7 and 8. Here again, the peak amplitude of 
the signal is not an increasing function of disbond and void 
size. In contrast, bond strength is an increasing function of 
the void and disbond sizes. 

tin the AE experiments, we expected to find a good cor- 
relation between AE parameters and bond strength for all 
types of specimens with different defects (steel-adhesive- 
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Fig. 14 Plots of the number of events along the lap joint for two typical specimens. 




Time, sec 

Fig. 15 Typical plot of AE energy vs. time for a) 
specimens with void or disbonds, b) perfect specimens. 


steel specimens, with 0 = 3.4). This is due to the nature of 
failure location and mechanism, and energy released during 
failure. Since for symmetric voids and disbonds in the cen- 
ter of the overlap, stress distribution is almost identical, the 
bond strength and energy released can be assumed to be the 
same for specimen with identical void or disbond size. 
Furthermore, all detected released energy is associated with 
the high stress areas, which are the two edges of the over- 
lap, which also control the bond strength. Figure 13 is a 
plot of a typical load-time curve of an adhesively bonded 
joint specimen along with its corresponding mean ringdown 
counts (number of threshold crossing) vs. time. There are 
three prominent peaks, the last of which corresponds to the 
time of failure. Having a lower amplitude, the first two 
peaks can be associated with crack initiation and prelimi- 
nary crack growth. The source of these prominent peaks, 
recorded in real-time, was found to be near the end of the lap 
joint where the failure was originated. Furthermore, most of 
the AE activities were confined to the edges of the overlap; 
see Fig. 14. 

Plots of energy vs. time compiled every four seconds 
during the period in which the specimens had reached their 
ultimate strength, exhibit two principal formats as shown 
in Fig. 15a. The specimens with a distinguished peak at 
failure, normally depicts specimens with artificial voids or 
disbonds. These specimens can be categorized as having 
failed in a brittle manner due to their short energy released 
duration. For these specimens, the energy released was rela- 
tively high at the peak compared to the perfect ones. 
Pollock and Lane (1968) showed that brittle materials 
exhibited high peak energy and short emission duration after 
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Fig. 16 Natural log of AE energy released at and im- 
mediately after the peak released energy vs. the square of 
maximum load (bond strength) 2 endured. 



Fig. 17 Natural log of the number of events at and after 
peak event (failure point) vs. the square of the maximum 
load (bond strength) 2 endured. 


the maximum load than ductile materials. This may be due 
to the fast fracture behavior in brittle materials and crack tip 
blunting (crack arrest) in ductile materials. For perfect spec- 
imens, fast fracture may have been prevented by having 
effective adhesive in the entire overlap. This is not true for 
specimens with void or disbond. Figure 15b depicts energy 
released for perfect specimens. Energy emission and the 
number of threshold crossings were sustained at a low level 
for some period after the time of the maximum peak of 
energy. This can be due to the crack or damage growth 
period in perfect specimens which is generally greater for 
the perfect specimens than specimens with a defect. 
Regression analyses of the total energy released and number 
of events at and after maximum peak energy show that 
these parameters can be related to the bond strength for all 
type of specimens; see Figs. 16 and 17. Analysis of the AE 
data further showed that the events with lower rise time 
were associated with the crack initiation while the one with 
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Fig. 18 Bond strength vs. minimum of ^weighted AUP) 
for aluminum-adhesive-aluminum specimens with various 
defects (mostly due to improper surface preparation and 
disbond). 

higher rise time was associated with the crack (damage) 
propagation. 

For the graphite epoxy composite-adhesive-graphite 
epoxy composite, 0 = 10. For this value of 6, it was 
expected that symmetric void or disbond up to 70% of the 
overlap length, to have a negligible effect on bond strength. 
Indeed there was little change in the bond strength for spec- 
imens with void up to 70% of the overlap length in the 
central region of the bonded joint (Nayeb-Hashemi and 
Rosseuos, 1990). The peak amplitude of the received signal 
in an SWF setup experiment (Fig. 9), was found not to be 
very sensitive to the bond strength and showed a trend simi- 
lar to the one exhibited by the dynamic analyses (Fig. 8). 
Other AU factors also failed to show any sensitivity to the 
bond strength. 

To address more effectively the effects of disbond, 
improper adherend surface preparation, over cured and under 
cured adhesive, and deficient pressure during manufacturing 
of these specimens on bond strength, we set up an immer- 
sion experiment (Nayeb-Hashemi and Rosseuos, 1993). As 
was described above, through transmission experiments uti- 
lizing a contact transducer resulted in little change in atten- 
uation data for specimens with disbonds. Furthermore, dis- 
bonds or voids may not result in any change in bond 
strength as shown for composites. Disbonds right at the 
edge and a disbond in the middle of the bond have com- 
pletely different effects on the bond strength. However, us- 
ing contact transducers may not be a sensitive means to 
evaluate any changes in attenuation. To incorporate the 
location into attenuation, peak amplitude, SWF and AUP 
measurements, a weighting function w(x) = exp(0x) was 
introduced. Here x is measured from the center of the over- 
lap. Using focus transducers with center frequency of 2 







MHz and adjusting transducer specimen distance, such that 
the focal point lies on the adhesive adherend interface, the 
peak amplitude, SWF and A UP at several locations of the 
verlap “re measured, using a pulse with center frequency 
2 Mi New ultrasonic parameters were introduced as 


Peak* = ^exp(0x i )Peak(x i ) 

(21) 

i-l 


SWF* = ^ exp(0Xj )S WF(x j ) 
i*l 

(22) 

AUP* = £ exp(0Xj ) AUP(x, ) 

(23) 


i*l 


These parameters were also evaluated by considering the 
weighting function to be equal to unity. These parameters 
were either evaluated for the entire overlap or the over lap 
was divided in half and these parameters were evaluated for 
each half. The bond strength was assumed to be controlled 
by the section with a lower specific parameter value. For a 
variety of aluminum-adhesive-aluminum specimens with 
many surface defects, minimum AUP produced an excel- 
lent correlation with bond strength (Fig. 18). However, the 
correlation of the same data with the bond strength without 
incorporating weighting function (minimum of or total of 
XAUP^) vs. bond strength) was not good. 

The results presented in this paper clearly show that the 
quest for an effective procedure for nondestructive evaluation 
of the adhesively bonded joint requires multidisciplinary 
approach. We believe that the immersion technique using 
focus transducers and incorporating the proper weighting 
function is the most effective method for bond quality inter- 
rogation. This procedure may be applied by using echo 
transducers. Further research has to be performed to substan- 
tiate the new parameters. 

8. Conclusions 

Defects such as voids, disbonds, improper adherends 
surface preparation and manufacturing procedures, may 
affect the bond strength of adhesive joints. In order to 
develop a proper acousto-ultrasonic method for bond quality 
interrogation, effects of defects such as void and disbond on 
(1) the stress distribution over the overlap, (2) the dynamic 
response, and (3) ultrasonic wave propagation must be 
understood. 

The stress distribution in a single lap joint with and 
without voids is analyzed. For identical adherends, the stress 
distribution depends on a nondimensionalized parameter, 0 . 
This parameter incorporates adherend and adhesive thick- 
nesses and properties. The analyses show that for joints 
with 6 > 10 the maximum shear stress is not affected by 
the presence of symmetric voids or disbond in the middle of 


overlap (void < 70% of overlap length). In contrast a finite 
element analysis shows that the dynamic response of a 
bonded joint is significantly affected by the presence of 
voids or disbonds in the overlap region. These results show 
that the measured ultrasonic parameters are affected by the 
defects. However, these measured values may or may not be 
related to the bond strength. 

Through transmission longitudinal and shear wave 
attenuation changes in steel-adhesive-steel specimens ( 0 * 
3.4) with and without defects (symmetric voids and dis- 
bonds in the middle of overlap region) were measured using 
pulse oscillators with the center frequency of 1 and 2 MHz. 
The experimental results showed that the specimens with a 
larger void exhibited higher attenuation. These specimens 
also had lower bond strength. In contrast, specimens with 
disbonds in the same range size as voids exhibited little 
variation in the measured attenuation change. However their 
strength was affected in a similar way as those with voids. 
Various other AU factors also failed to show good sensitivi- 
ties to the bond strength. 

Since the location of defects has an important effect on 
the shear stress distribution and subsequent bond strength, a 
weighting factor w(xj) = exp( 0 Xj) was suggested. Using an 
immersion system with focused transducers many AU fac- 
tors were evaluated. New ultrasonic parameters based on 
measured values and weighting function were proposed. It 
was shown that some of the new parameters produced a 
good sensitivity for bond quality prediction. However the 
sensitivity of the same parameters without weighting func- 
tion was somewhat poor. 

AE activities during the tensile tests were monitored 
and several time-domain variables were gathered for steel- 
adhesive-steel specimens. Two different formats of energy 
vs. time were observed. The specimens with a perfect bond 
showed an abrupt energy release at failure, followed by a 
short period of continuous AE activities. However, speci- 
mens with voids or disbonds also showed abrupt energy 
release, but followed with no additional activities. This 
behavior was related to the crack initiation and propagation. 
AE parameters were found to be a means for bond quality 
prediction. 
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The only viable method for joining plastic tubes and composite shafts is by bonding them adhesively. 
These structures are often subjected to complex cyclic loadings. Failure of these tubular joints not only 
depends on the applied loads, but also depends on the tube geometry, material properties of adhesive 
and tubes, and defects in the joint. 

The shear stress distribution in the tubular joints is obtained for joints under axial and torsional 
loadings using the shear lag model. Under axial loading the adhesive is assumed to carry only shear 
stress and adherend to carry only axial load. However, the model considers the variation of the shear 
stress across the adhesive thickness. The effect of a void on the maximum shear stress is obtained. 

A nondimensional 0 a parameter is defined and it is shown that the shear stress distribution not only 
depends on the value of 0 a but it also depends on cross sectional geometry of the tubes. For tubes with 
equal cross sectional area, the shear stress distribution along the bonded area is almost symmetric. For 
tubular joints with 0 a equal or greater than 6.7, a centrally symmetrical annular void with a size of at 
least 50% of the overlap length has little effect on the maximum shear stress and thus the failure load. 

The shear stress under torsional loading is obtained by assuming the adhesive to shear in the 
circumferential direction only and neglecting its other deformations. The tubes are assumed to shear in 
the axial direction. The analysis considers the variation of the shear stress across the adhesive thickness. 
As in the case of axial loading, a new nondimensional parameter, 0 t , for tubes under torsion is defined. 
The results show that the shear stress in the bonded area not only depends on the 0 t value, but also 
depends on the polar moment of inertia, J\ and of the tubes. The effect of annular voids on the shear 
stress distribution is evaluated for different void sizes and 0 t values. 

The failure locus of adhesively bonded tubular specimens under axial, torsional and combined axial 
and torsional loadings is obtained. Based on these results a damage model for the tubular joints under 
combined axial and torsional cyclic loading is proposed. It is shown that this model can predict the 
fatigue life of the tubular joints reasonably well. © 1997 Elsevier Science Limited. All rights reserved. 


(Keywords: tubular joints; composites; destructive testing; stress analysis; fracture mechanics; destructive testing by lan- 
shear joints) 6 J ^ 


INTRODUCTION 

Composite and polymeric tubes are becoming more 
popular in structural applications and transmission 
shafts. These tubes are often connected to each other 
with an adhesive. The bond strength and its 
degradation during service depend on the mechanical 
properties of the tubes and adhesive, the geometry of 
the tubes, defects present in the bonded area, and 
loading conditions. There are many theoretical analyses 
addressing the stress distribution in the tubular joints 
under axial and torsional loadings. However, few of 
these theoretical works have considered the effects of 
defects, such as a void in the adhesive, on the shear 
stress distribution. Furthermore to the best of our 
knowledge, there is no literature on the multiaxial 
fatigue life estimation of tubular joints under combined 
axial and torsional loadings. 


Lubkin and Reissner have analyzed the stress 
distribution in tubular lap joints under axial load and 
gave solutions for both the shear stresses, r zr , in the 
adhesive layer and the normal stress, <x n , across the 
thickness of the adhesive layer, due to adherend 

bending 1 . In the analyses, the two adherends were 

assumed to be thin, and hence the thin shell 

approximation was employed. The adhesive layer was 
assumed to be thin and much more flexible than the 
adherend, and was modeled as a series of infinitesimal 
tensile and shear springs. The work of the stresses, r zr 
and <r rr , in the adherend is assumed to be negligible 
compared to the work of these stresses in the adhesive. 
The effect of defects in the analyses was also neglected. 

The stresses in adhesively bonded tubular lap joints 
subjected to axial and torsional loads have been 

analyzed using axisymmetric quadratic isoparametric 
finite elements by Adams and Peppiat 2 . In the axial 
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where i a is the shear stress in the adhesive. Using the 
adhesive constitutive equation, equation (3) can be 
written in terms of adhesive displacement w a as 


^ Wa , du a 

r d7“ U 


(4) 


where <J a is the shear modulus of adhesive. Assuming 
the displacements of adherend 2 and 1 are u 2 and iq at 
r = R 3 and r — R 2 > respectively, equation (4) can be 
solved to obtain the adhesive displacement field across 
its thickness. This can be presented as 


Wa = W 2 + 


U 2 “ Wj 


In R 3 — In R 2 


(Inr - \nR 3 ). 


(5) 


Using equation (5) the shear stress at the inner and 
outer surfaces of adherend 2 and 1 can be found from, 


Orx) 2 = G. 


d«a 

"dT 


dw a 


= U a 


(T«)i = 




r — ^2 


W 2 - M, 


[In R 3 — In R 2 \ R 


1 


IT ( 6 ) 


w 2 - w, 


In R 3 — In R 2 


1 

Ri 


(7) 


Substituting (6) and (7) into (1) and (2) results in 
d 2 « 2 


and 


W-Ri)^-r-2G a 




U 2 - M, 


In R 3 — In R 2 


u 2 ~ U\ 


\nR 3 - ln/L 


-0 ( 8 ) 


-0 (9) 


Letting x = where L is the length of the overlap, 
equations (8) and (9) can be written as 


EAtf - R*) d 2 m 2 


U 


d? 


and 


£,(**-**) d 2 M 2 

L 2 d { 2 


— 2G, 


-2G a 


W 2 - U\ 


In R 3 — ln/? 2 


w 2 - U 1 


[In /? 3 - In /? 2 J 


Equations (10) and (11) are valid for the bonded 
region of the joint when an adhesive layer is present. 
For the region without an adhesive layer, the 
equilibrium equations are 


and 


d 2 w 2 

W 

d 2 M| 

w 


— o 


= 0 


( 12 ) 


(13) 


The solution of equations (10) — (1 3) gives displacements 
in different region of the overlap. The displacement 
fields in the region 1 where the adhesive layer is 
present, are 


u 2 (0 = C 3 e^ + C 4 e- fl * <: - g(C| ^ + C2) (14) 


«.«)= 


i 


EiiRl-R]) 


[QZ + C 2 -E 2 (Rl-R 2 3 )u 2 ] 


and for the region (2) with a void are 

(15) 

MO = C 5 { + C 6 

(16) 

«,(£) = C 7 £ + C 8 

(17) 

and for region (3) where again an 
present are 

adhesive layer is 


«2(0 = C 9 e e ' ( + C l0 ^ e ‘ ( - ^ (C„<f + C l2 ) (18) 


and 


wi (<D — [C^ + ^12 “ £2(^4 ~ ^3)^2] 


1 


£,(/? 2 -/??) 


(19) 


where 0 a is a nondimensional parameter 
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, ! ^(^-^ 
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d continuity conditions in 
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the various 
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and 



d 2 u 2 

dx 2 


- 2nG a 


u 2 - u i 

J 1_ 

R 2 R 3 


= 0 


(46) 


Equations (45) and (46) can be normalized by letting 
x ~ ZL and then solving for u x and u 2 in terms of 
This results in 


u 2 (£) = C 3 e e,i + C 4 e~ e,i + C 2 ) (47) 

and 




r 2 l 2 


+ c 2 — 


J 2 G 2 

L 2 R 3 


m 2 


(48) 


where here again is a nondimensional parameter 
for the case of torsional loading, which combines 
adherend and adhesive mechanical properties and 
geometries, 



R 2 RlL 2 


(49) 


and B is a constant depending on adhesive and 
adherends mechanical properties and geometries 




-2 nG a 


1 R\R\l? 

Ri - Ri {J\G,)J 2 G 2 


(50) 


Equations (47) and (48) are valid in region of a tubular 
joint when adhesive layer is present. For a region with 
an annular void the rotational displacements and w 2 
are 


k 2 = C 7 £ + C 8 (51) 

Ml = C 5 £ + c 5 (52) 

A similar procedure as that developed for axial loading 
can be used to obtain rotational displacements at 
various regions of the overlap. The boundary 
conditions are also similar to those for the axial case 
except that at £ — 1 


du 2 TLR , 

The shear stresses (z r0 ) w and (r r0 ) 2} can then be 
obtained by substituting for u x and u 2 from equations 
(47) and (48) into equations (40) and (41). These 
stresses are evaluated for various 9 and J x and J 2 
values. The normalized shear stresses are defined as 
the ratio of the actual stress to the average stress, 
giving 


SCF = 


( T r0)2 i fafl) li 


( T av)2i ( T av)l0 

where (i av ) 2i and (t av ) J0 are 

( T av) 2 i — and (t -)-o = OTZ 


(54) 


(55) 


EXPERIMENTAL INVESTIGATION 

Figure 3 shows the specimen geometries used to 
investigate the failure locus and multiaxial fatigue 
properties of tubular bonded joints. Steel adherends 
were bonded using an adhesive with 50% epoxy and 
50% hardener supplied by the Shell Company (Epoxy 
828 and Hardener V40). The mixture was degassed for 
3 min prior to its usage. For this mixture of the 
adhesive, the shear modulus of the adhesive was 
reported to be 791 MPa 18 . The adhesive was cured at 
100°C for lh. The surfaces of both adherends were 
sand blasted and cleaned ultrasonically in acetone. This 
resulted in highly repeatable experimental data. A 
fixture was used to ensure alignment of adherends 
prior to adhesive curing. A rubber gasket was used 
between male and female adherends to retain adhesive 
during its curing and avoiding its leakage, Figure 3. 

For the failure locus experiments, specimens were 
installed in an Instron tension/torsion machine (Model 
1322) and pulled or twisted to failure at the rate of 
70 N or (2 Nm)/s. Increasing or decreasing the loading 
rate two or three times did not result in any change in 
the failure loads. An analog/digital convertor was used 
to measure Ioad/displacement. The data were stored in 
a PC for future analysis. For fatigue experiments, the 
specimens were subjected to sinusoidal loading with 
zero mean load at the frequency of 1 Hz. The 
experiments were conducted under constant 
displacement amplitudes. The failure was defined when 
either the axial load or the torque dropped by 10%. 
Most of the failure was initiated at the inner surface 
where stresses were maximum. 



139.7- 


i— - 


-139.7- 


L 


24.0 


50.8 


-200.7- 


391.2 



All Dimensions are in nn 

Figure 3 Schematic diagram of the steel adhesive- steel tubular joint 
used in the experimental investigation 
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Normalized Distance, X/L 

Figure 8 Normalized shear stress, 7^/^, distributions in a tubular 
bonded joint under torsion, with and without an annular void. The 
void is located at the edge where the torque is applied on the 
adherend with a smaller J value 



Normalized Distance, X/L 


Figure 9 Normalized shear stress, T r0 /T avc , distributions in a tubular 
bonded joint under torsion, with and without an annular void. The 
void is located at the edge where the torque is applied on the 
adherend with a larger J value 



Normalized Distance, X/L 

Figure 10 Normalized shear stress, T^/T ave , distributions in a tubular 
bonded joint under torsion, with a central annular void of different 
sizes 



Normalized Distance, X/L 

Figure 11 Normalized shear stress, T rjc /T ave , distributions in a tubular 
bonded joint under tension, with a central annular void of different 
sizes 



Figure 12 Normalized shear stress, (x r;c ) ma!l /T avc , versus normalized 
void size in tubular joints under tension for several values of 9 



Shear Stress x /(% ) 

rx' v n^max 

Figure 13 Failure locus of a steel-adhesive-steel tubular joint under 
in phase axial/torsional loading 
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Theoretical and Experimental 
Evaluation of the Bond Strength 
Under Peeling Loads 

Reliable applications of adhesively bonded joints require understanding of the stress 
distribution along the bond-line and the stresses that are responsible for the joint 
failure. To properly evaluate factors affecting peel strength, effects of defects such 
as voids on the stress distribution in the overlap region must be understood. In this 
work ; the peel stress distribution in a single lap joint is derived using a strength of 
materials approach. The bonded joint is modeled as Euler-Bemoulli beams ; bonded 
together with an adhesive, which is modeled as an elastic foundation which can resist 
both peel and shear stresses. It is found that for certain adhesive and adherend 
geometrtes and properties, a central void with the size up to 50 percent of the overlap 
length has negligible effect on the peak peel and shear stresses. To verify the solutions 
obtained from the model, the problem is solved again by using the finite element 
method and by treating the adherends and the adhesive as elastic materials. It is 
found that the model used in the analysis not only predicts the correct trend for the 
peel stress distribution but also gives rather surprisingly close results to that of the 
finite element analysis. It is also found that both shear and peel stresses can be 
responsible for the joint performance and when a void is introduced, both of these 
stresses can contribute to the joint failure as the void size increases. Acoustic emission 
activities of aluminum-adhesive-aluminum specimens with different void sizes were 
monitored. The AE ringdown counts and energy were very sensitive and decreased 
significantly with the void size. It was observed that the AE events were shifting 
towards the edge of the overlap where the maximum peeling and shearing stresses 
were occurring as the void size increased. 


Introduction 

Adhesive bonding has been used extensively in the aerospace 
and other high-technoiogy industries and has a great potential 
for applications in other areas of manufacturing. It is attractive 
because it distributes stress over the entire bond area and elimi- 
nates the stress concentrations which can occur with mechanical 
fasteners. Bonded joints have potential advantages of strength- 
to-weight ratio, design flexibility and ease of fabrication. 

A Fokker Aerospace report [1] showed that many of the past 
failures of bonded structures involved poor design. Most often 
the poor design could be attributed to inadequate understanding 
of the adhesion failure mechanisms. 

There are three dominant modes of mechanical failure in 
bonded joints: 1) adherend failure (including delamination if 
one of the adherends is a laminated composite), 2) interfacial 
failure between adhesive and an adherend and 3) cohesive fail- 
ure within the adhesive. When an adhesively bonded joint is 
subjected to an out-of plane loading, the tearing of the adhesive 
which occurs is called peeling. Because this type of failure 
can be produced by normal loads which are relatively small 
compared to the shear loads which structural adhesives are capa- 
ble of withstanding [2], peel strength of adhesively bonded 
joints is a property to be considered. 

In its simplest form, an adhesively bonded structure consists 
of three components of different mechanical properties, namely 
the adhesive and the two adherends. Under most operating loads 
and environmental conditions, the adherends behave in a lin- 
early elastic manner, however, the adhesive may exhibit visco- 
elastic or nonlinear behavior. The exact analytical solution of 
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the problem regarding the stress distribution in the bonded area 
is complicated. The existing analytical studies are, therefore, 
based on certain simplifying assumptions with regard to the 
modeling of the adhesive and the adherends. The adherends are 
usually modeled as an isotropic or orthotropic membrane [o], 
plate [4], or elastic continuum [5]. The primary physical con- 
sideration used in the selection of a particular model is generally 
the adhesive-to-adherend and adherend- to- adherend thickness 
ratios and the ratio of the adherend-to-adherend thickness to its 
lateral dimensions. 

Erdogan et al. [6] have analyzed a general plane strain prob- 
lem of adhesively bonded structures which consist of two differ- 
ent orthotropic adhereuds. The thickness of the adhesive was 
assumed small compaied to the thickness of the adherends 
which, in turn, are small compared to the length of the joint. 
The transverse shear stress effects in the adherends and the in- 
plane normal strain in the adhesive are taken into account. The 
solution is obtained by assuming linear stress-strain relations 
for the adhesive. The peak values of the shear as well as that 
of the normal stress in the adhesive are found to be at the 
edges of the overlap region. In another study [7], an adhesively 
bonded lap joint is analyzed by treating the adherends as linear 
elastic plates and the adhesive as linearly viscoelastic solid. The 
stress distribution in the adhesive layer is calculated for three 
different external loads, membrane loading, bending, and trans- 
verse shear loading. The results indicate that the peak value of 
the normal stress in the adhesive is not consistently higher than 
that of the corresponding shear stress, and its distribution decays 
slower than that of shear stress from the edge of the overlap. 
In [8], the governing equations for a step lap joint with a void 
are established using a modified shear lag model, where the 
adhesive can have extensional as well as shear deformations. 
The model considers a quadratic axial deformation across the 
adhesive thickness. This model was used in order to accommo- 
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date thermal mismatch between the adhesive and the adherend, 
where the adhesive can expand due to temperature changes. It 
was shown that relatively large void sizes have little effect on 
peak shear stress for sufficiently large values of one of the 
defined parameters, and the value of this parameter can be con- 
trolled by either changing the geometry of the bonded joint or 
its materials. 

There are a substantial number of papers dealing with bond 
strength prediction using nondestructive evaluation techniques. 
Acoustic emission is a nondestructive evaluation technique that 
has been known to have a great success in predicting structural 
integrity of components. 

Acoustic emission (AE) is the generation, propagation and 
detection of transient stress waves in materials as they undergo 
deformation or fracture. These waves propagate to the surface 
of the structure where they may be detected by an ultrasonic 
transducer. The output of the transducer is processed and the 
resulting signals are interpreted as the ‘‘AE” signals. Acoustic 
emission signals can take many forms depending upon the mate- 
rial and the failure mechanism in the material. AE signals from 
defects in composites and geological materials generally contain 
information at low frequencies 0.5-100 KHz, since attenuation 
is relatively high due to the complexity of these materials. Sig- 
nals of significance in metals and brittle materials contain infor- 
mation between 100 KHz and 2 MHz. In this range, a good 
compromise is found for most testing applications because am- 
bient noise is low [9] . 

A number of techniques are employed to isolate and validate 
signals from noise in the time domain. Times of arrival can be 
used to permit geometric elimination of obvious noise through 
gating, and acceptance of only those signals which arrive from 
a particular region of the structure. In time domain, acoustic 
emission data include: ringdown counts, rms voltage, number 
of events, energy rate, rise time, event duration, amplitude distri- 
bution, and numerous others. 

Williams and Lee [10] in a comprehensive study of acoustic 
emission in fiber composites materials and structures, found 
that AE is capable of detecting potential failure sources and 
defects in complex structures independent of the location, type 
and orientation of the flaws. They also stated that in tests involv- 
ing AE, both the geometry and materials may play an important 
role in the AE results. Accurate and reproducible results require 
calibrations of specimens or structures used. Hashemi et al. [9] 
monitored acoustic emission activities during tensile tests of 
steel single-lap joint specimens. Several time domain parame- 
ters were gathered and analyzed. They observed two different 
formats of energy versus time. The specimens with a perfect 
bond showed an abrupt energy release at failure, followed by 
a short period of continuous acoustic emission activities. How- 
ever, specimens with voids or disbonds showed abrupt energy 
release, but followed with no additional activities. They reported 
that this behavior was related to the crack initiation and propaga- 
tion and that acoustic emission parameters were found to be a 
means for bond quality prediction. Williams et al. [11] moni- 
tored acoustic emission in adhesively bonded automotive glass 
fiber composite single-lap joints. Their AE data obtained during 
monotonic loading was analyzed using the AE load delay con- 


cept. The concept of the “AE stress delay” as defined by Wil- 
liams et al. [12-14] is used to characterize the AE-strength 
behavior of the specimens. The AE stress delay is defined as 
the stress required to produce a specified level of cumulative 
AE events. They found a linear relation between the AE load 
delay and the fracture strength of the bonded joints, independent 
of the flaw type in the specimens. Also, the AE load delay is 
found to be useful in distinguishing the delamination, interfacial 
and cohesive separation modes. The purpose of this research is 
to understand the effect of defects on the peel stress distribution 
and nondestructive evaluation of peel strength of bonded joints 
by acoustic emission and acousto-ultrasonic technique. A 
strength of materials approach and the finite element technique 
are used to find the stress distribution in single lap joints sub- 
jected to a direct peeling load. 

Theoretical Investigation 

The stress distribution in a single lap joint is derived using 
a strength of materials approach by assuming that the adherends 
behave as Euler-Bemoulli beams. The adherends in the bonded 
region are assumed to be supported by a continuous elastic 
foundation. When the adherends are deflected, the intensity of 
the distributed reaction at every point of the overlap is taken to 
be proportional to the difference in the deflections of the ad- 
herends at that point. Under such conditions, the reaction per 
unit length of the adherend can be presented as ky, where y is 
the relative deflection of the adherends at point x, k is the 
stiffness of the foundation and its value is obtained from k = 
(E a w/r) where E a is the elastic modulus of adhesive layer, w 
is the width of the overlap, / is the adhesive thickness. In previ- 
ous work [15], the adhesive layer was modeled as a combina- 
tion of linear and torsional springs. The torsional spring constant 
of the adhesive was obtained by modeling the adhesive as a 
cantilever beam. The torsional spring constant k n depended on 
a scalar factor ranging from 0 to 1. In order to improve that 
model, the effect of the shear stress at the adhesive-adherend 
interface is taken into consideration without any assumption 
regarding its value. Figure 1 shows a single-lap joint geometry 
and dimensions. For a perfectly bonded joint with no void, the 
structure is divided into three regions, where the lateral deflec- 
tion of adherend in region 1 and 2 is taken as y x , and y 2 and 
that of adherend 2 in region 2 and 3 are taken as y 3 , and y 4 . 
Let V, A/, and r be the transverse shear load, moment and shear 
stress in the adhesive respectively. From the free body diagram 
in Fig. 2, the equilibrium equations in the region 2 can be 
written as: 


dV 

— + k(y 2 - yj) = 0 
dx 


V + 


rwh I 

~r 



Substituting for V from Eq. ( I ) into Eq. (2) results 

+ _ 0 
dx 2 dx dx 1 


( 1 ) 

( 2 ) 


(3) 
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Fig. 2 Free body diagram of adherend 1 and adhesive 
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Region 3, adherend #2 


<i J v. 

(Elh- i = 0 
dx 


The bending moment and deflection are realted as 

d 2 v-> 

M = (El), -pr 
dx~ 

Substituting for M from Eq. (4) into Eq. (3) yields 
d^V-y wh\ dT 

(E/)i — — + k{yi — yj) = 0 

dx 2 dx 


Boundary 1 conditions: 
• Ac x = 0 


d } Vj 

E/ ' = F 
(Elh 4^ = 0 

dx~ 


At * = /, 


The Shear strain y ad and shear stress r at the adherend-adhesive 
interface can be determined by assuming that the shear strain 
in the adhesive is proportional to the relative slope of adherend 
1 with respect to adherend 2. This can be written as: 


h\ dy 2 h 2 dy 2 

7 "7" ~ 77T 


The shear stress r can then be obtained from: 


r _ G h\ dy x hi dy 2 
T 7<ui t 2 dx 2 dx 


Where G is the shear modulus of the adhesive, A t , and h 2 are 
the thicknesses of adherend l and 2 respectively and (dy 2 /dx) 
and (dy 2 /dx) are the slopes of adherend 1 and 2. Substituting 
for r from Eq. (7) into Eq. (5) results in a general form of a 
beam on an elastic foundation with a combined tensile and shear 
resistances. 


d*y 2 

II- 
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dx 4 
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The modulus of bending rigidity of the adherends is denoted 
by (£7) ( and ( £7) 2 . The governing deflection equations for the 
entire bonded joint with one end clamped and the other end 
subjected to a peeling force F can then be written as: 


Region 1, adherend #1 


(Elh = 0 

dx 


Region 2, adherend #1 
d A y~> 

(E/), + k(y 2 - y 3 ) 


d 1 v. 

(Elh ~fr = 0 

dx' 
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Fig. 3 Schematic diagram of the deformation field in the bonded joint 
where the shear strain was derived from 
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Fig. 4 Schematic diagram of the acoustic emission data acquisition system 
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Similar governing equations and boundary conditions can be 
written for the bonded joint with a void, by dividing the overlap 
into three regions. This analysis was performed and the effect 
of voids on the peel and shear stresses was derived. 

In order to verify the results obtained for the peel and shear 
stresses from the model, these stresses were again obtained by 
finite element analyses, using ADINA finite element code. The 
bonded joint was modeled as 2 dim. structure, using 8 node 
isoparametric elements. The analysis was done for the same 
materials and geometry used in the experimental investigation. 
Voids were created by giving the modulus of the adhesive for 
elements in the void region a very small value. Adherends were 
6061-T6 aluminum with elastic modulus of E - 69 GPa. The 


analyses were performed for joints with adherends of equal or 
different thicknesses. The adhesive layer was Hysol EA9689 
adhesive with E = 2.2 GPa and a thickness of .13 mm. The 
overlap was 25.4 mm x 25.4 mm and symmetric central voids 
were introduced in the overlap. Effects of the mesh size in 
the overlap region were also investigated. Based on the results 
optimum mesh size was selected in the overlap region. 

Experimental Investigations 

Adhesively bonded joint specimens were prepared using alu- 
minum 6061-T6 as adherends. Adherends were joined together 
in a single lap joint configuration using Hysol EA9689 epoxy 
film of .13 mm thickness. The adherend surfaces were sand- 
blasted and cleaned with acetone prior to bonding. The overlap 
area was 25.4 mm x 25.4 mm and upper and lower adherend 
thicknesses were 3.17 mm and 12.6 mm, respectively. Voids 
were introduced in the overlap by cutting and removing the 
adhesive film from the desired area of the boned joint and 
spraying a mold release agent over the void area before bonding. 
The specimens were cured at 170°C for 1 h, and dummy tabs 
and shims were used during the curing process to ensure that 
the adhesive layer in the specimens had the same thickness. 
Strain gages were also mounted on the surface of each specimen 
in order to monitor strain changes up to the specimen failure. 
Identical specimens were manufactured in order to ensure re- 
peatability of the obtained data. An MTS model 810 servo- 
hydraulic testing system was used for the test. 

The data for acoustic emission were monitored and collected 
using AET 5500 system. Two AE transducers, each with a 
center frequency of 175 KHz and nominal diameter of 22 mm, 
were used to detect events emitted from the specimens at a 
distance of 76 mm (center to center of transducers 98 mm) 


Flexible Adherends Flexible Adherends 




Fig. 5 Comparison of the shear stress distribution obtained from the Fig. 6 Comparison of the peel stress distribution using theoretical 
theoretical model and finite element model and finite element 
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Fig. 7 Normalized variation of the maximum peel stress as a function 
of the void size 


apart. Each sensor was coupled with the workpiece through a 
B-type panametrics couplant and held in position with the aid 
of 4 stiff springs as shown in Fig 4. Data collection were focused 
on the lap joint portion only and all the other AE activities were 
thus discarded. AE event locations were established using the 
difference in the arrival times of the signals received by the 
two sensors. ASTM standard of pencil lead fracture method 
[24] was used in order to calibrate the system prior to the actual 
run. Calibration eliminated the need to obtain the velocity of 
emitted waves. It is assumed that lead pencil fracture simulate 
stress waves generated by the actual failure and first arival wave 
is the longitudinal wave. The detected signals were preamplfied 
using two 60 db preamplifiers and a threshold of 0.5 V was 
chosen to eliminate the background noise. The entire wave 
forms and important features of the AE wave forms were re- 
corded using a P.C. After completion of the tests, data were 
analyzed and processed using an AET 5500 software (Ba- 
win™) and various AE parameters were extracted. 

Theoretical Results and Effect of Voids on Stress Dis- 
tribution 

Figures 5 and 6 show the shear and peel stress distributions 
obtained from the theoretical and finite element analyses. The 
results show good agreement for shear and excellent agreement 


c 

o 



Fig. 8 Normalized variation of the maximum shear stress as a function 
of the void size 



for peel stresses between the theoretical and finite element re- 
sults. In our previous work, the peel and shear stress distribution 
in the bonded joint were obtained by modeling the adhesive 
layer as a combination ot linear and torsional springs. The tor- 
sional spring constant was determined from a cantilever beam 
model which resulted in a spring constant with a scalar factor 
[15]. Although the model was capable of determining the peel 
stress, it resulted in a poor determination of the shear stress at 
the adhesive/adherend interface. Furthermore, the scalar factor 
for the torsional spring constant was determined from the solu- 
tion obtained from the finite element results. The current analy- 
ses overcomes these deficiencies and provides excellent results 
for both peel and shear stresses without any assumptions and 
the need for a scalar factor. The model was further verified for 
different bonded joint dimensions and again good agreement 
was found between peel and shear stresses and those obtained 
from the finite element analyses. Results show that when both 
adherends are flexible and have the same thickness, both peel 
and shear can contribute to the failure of the joint when a void 
is introduced. But, when one adherend is flexible and the other 
is rigid, peel stress become dominant and contribute mainly to 
the failure. The effect of the void in the adhesive was also of 
interest in this investigation. It was found that for certain adhe- 
sive and adherend geometries and properties, a central symmet- 
ric void had a negligible effect on the peak peel and shear 
stresses. Figures 7 and 8 indicate the variations of the maximum 
normalized (normalized with respect to the maximum stresses 
for the case of no void) peel and shear stresses for a joint with 
a central void size up to 70 percent of the overlap length. Figure 
9 show the effect of a central void on the peel stress distribution 
in the joint. These results are similar to those obtained for single 
lap joints under tension and tubular joints under tension/torsion. 



Fig. 10 Variation of the damage parameter with the void size usinq finite 
element results 
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For these joints it was also found that the introduction of central 
symmetric void may have little effect on the bond failure [9, 
25]. These results could not justify experimental data since 
significant changes in the failure loads were observed with the 
introduction of central voids in the bonded joint. 

In order to obtain a failure criterion for the bonded joint, a 
damage parameter D - aa 2 + 0r 2 , where a and r are peel 
and shear stress at the tensile end joint was defined. The bonded 
joint is assumed to fail when D = 1. This concept is analogous 
to the failure criteria such as yielding. A material is considered 
to fail if D > 1. Similar models have been developed for fatigue 
life prediction of components under multiaxial loadings [25]. 
The values for a and 0 were obtained from the failure load on 
two specimens (a perfect specimen with no void and a specimen 
with 50% central void) and using the stresses obtained from 
the finite element results. After establishing the values of a and 
0 the damage parameter D were found for different void sizes. 
Figure 10 shows that the damage parameter monotonically in- 
creases as the void size increases. It may be thus concluded 
that this is a better criterion to relate NDE parameters to the 
bond strength than the individual peel or shear stresses. Indeed 
our obtained AE data support this conclusion. 

Experimental Results 

Four specimens, one without a void and three with 6.3 mm, 
12.6 mm and 20.2 mm wide voids were subjected to an increas- 
ing peeling load up to failure. At the same time, strains readings 
were recorded and acoustic emission activities were monitored. 
All detected events are associated with the high stress areas, 
which is the tensile edge of the joint which also controls the 
bond strength. Figures 11, 12 show that the Cumulative AE 



Fig. 12 Variation of AE cumulative energy with the void size 
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Fig, 13 Variation of peak amplitude with the void size 


ringdown counts decreases with the void size, also that the AE 
energy exhibit the same behavior. Thus, these two parameters 
appear to be a good indication of the damage brought to the 
joint by the presence of the void in a set of specimens that 
otherwise are identical. The average of events peak amplitude 
is also sensitive to the presence of the void and increases with 
the void size as shown in Fig. 13. It is also found that the 
AE events are shifting towards the tensile end joint where the 
maximum peeling and shearing stresses are occurring as the 
void size increases as shown in Fig. 14. Further work is in 
progress to evaluate acousto-ultrasonic parameters for various 
bonded joints and relating them to the peel strength of the joint. 

Conclusion 

The stress distribution in a single lap joint under direct out- 
of plane loading with and without voids is analyzed by modeling 
the adherends as Euler-Bemoulli beams on an elastic founda- 
tion. The results show that the maximum peel and shear distribu- 
tions are confined to the edge ends of the overlap. The results 
also show that for the type of adherend and adhesive analyzed, 
a void up to 50 percent of the overlap length has a negligible 
effect on the peak peel and shear stresses. The stress distribution 
in the joint is compared to the results obtained from the finite 
element analyses. The theoretical model and finite element re- 
sults are in good agreement. A damage parameter was defined 
and shows promise as a criterion to relate NDE parameters to 
the bond strength. 

Acoustic emission activities were also monitored and several 
parameters were extracted. Cumulative AE ringdown counts 
and energy were found to be sensitive and decreased signifi- 
cantly with the defect size. 



Fig. 14 Distribution of events by location in specimens with different 
void sizes 
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Dr. Alex Vary December 2, 1994 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

Dear Alex: 

Enclosed you will find our recent publications made possible by your continuous support. 
Regarding our current efforts, I will outline some of our progress for your information. Because our 
work is still progressing, I am requesting a no cost three month extension for our project. 

Nondestructive Evaluation of Adhesively Bonded Joints 

Our efforts have continued toward understanding factors affecting the tensile strength of single 
lap joints in tension as well as peel strength of bonded joints and their effects on the acousto ultra- 
sonic parameters. The peel stress distribution in the bonded joint was obtained by modeling the single 
lap joint as a beam on an elastic foundation. Fig. 1. As with the case of the shear stress distribution in 
the single lap joint, the maximum peel stress was confined again to an area close to the edge of the lap 
joint. Central voids again had little effect on the maximum peel stress. Thus one may conclude again 
that the quality of the bond under peeling is controlled by the material close to the edge of the over lap, 
Fig.2. We have also studied wave propagation in the lap joint using finite element techniques, Fig. 3. 
This work was undertaken to address comments received from some reviewers, since the previous 
analysis was based on the steady state response of the lap joint rather than the transient response. 

The stress wave finite element model resembled the acousto ultrasonic set-up. A broad band 
pulse with the center frequency of 2 MHz was introduced into the specimen and the received wave 
was analyzed. Fig. 4. Here as in the case of the vibration analysis reported earlier, the peak amplitude 
of the received signal decreased with increasing defect size. Fig. 5. However, as reported here and 
earlier, the peel and shear strength did not exhibit a similar trend. Thus, our previous conclusions hold, 
and using a weighting fiinction in conjunction with the measured acousto-ultrasonic parameters 
appears to be justified. 

We have also been involved in nondestructive evaluation of tubular bonded joints subjected to 
tension and torsion. The preliminary work is concentrated on nondestructive evaluation of adhesively 
bonded tubular aluminum tubes. Future work will consider tubular composite joints. The failure locus 
of aluminum joints under tension and torsion is found. Fig. 6. The effects of defects in the adhesive 
region on the tension and torsion strength of the joint and its relation with the acousto ultrasonic 

334 Snell Engineering Center 
N'orrheastem University 
Boston, Massachusetts 02115 
617-373-2982 (office) 617-373-2921 (fax) 



properties are the focus of our investigation. The degradation of the bond during cyclic 
tension/torsion fatigue is another area of investigation. 

Nondestructive Evaluation of Composite Materials Subjected to a Localized Heat 
Damage (Burning) 

Unidirectional fiber glass epoxy composite panels with a thickness of 0.86 mm 
were manufactured and 12.7 mm strip composite specimens were cut from it. The mid 
point of the specimens was then subjected to 400 C by pressing a 6.35 copper cylinder at 
400 C against the composite surface. The contact time was varied from 5 sec. to 9 min. 
Since the glass transition temperature of the E glass fiber in the composite was much 
higher than 400 C, most of the damage was associated with the matrix decomposition. 
The extent of the matrix decomposition depended on the exposure time. However, it 
seemed that after 2 min. exposure, the entire 'matrix in the mid section was decomposed 
and further exposure caused the decomposition of the matrix along the specimen length. 
The specimens were subsequently interrogated using an acousto ultrasonic set-up. A 
broad band pulse with the center frequency of 2 MHz was used as the transmitting signal. 
A resonant Panametric transducer with the center frequency of 2 MHz was used as the 
transmitting transducer. The receiving transducer was a Fc-500 AET transducer with a 
flat sensitivity in the frequency range of 100 kHz-2 MHz. Stress wave factor, SWF, and 
acousto ultrasonic parameter, AUP, were measured for these specimens. The specimens 
were then pulled to their failure, and acoustic emission activities of these specimens were 
measured. The results showed that there is no significant strength reduction beyond 2 
minutes exposure at 400 C. This may be explained by the fact that the strength of the 
composites is controlled by the fiber strength. AUP of these specimens also showed a 
similar trend as those strengths versus the exposure time, figures 7, and 8. This may 
again be justified since the major part by of the stress wave energy is transmitted along the 
fiber. This phenomena is also reported in the literature. Specimens without heat exposure 
and those with exposure more than 1 minute showed different failure pattern and different 
acoustic emission activities. Figures 9 and 10 show events per location of perfect 
specimen and specimen with one minute heat exposure. The results' show that most 
activities of the specimens are concentrated in the mid section of the specimens, while 
events from specimens with no heat exposure are distributed along the entire gage 
section of the specimen. Acoustic emission activities of specimens with two or more 
minutes of heat exposure were also relatively unchanged. We are currently analyzing our 
data and I will inform you about our results. This is an important area of research, since 
many composites may be subjected to localized burning due to malfunctions of electron- 
ic components attached to composites or heat from other sources. 

A preliminary finite element stress analysis of this problem has also been 
performed, by conducting the static analysis of a long composite plate with a damaged 
region. The damaged region, where the matrix has deteriorated and the fibers are still in 
tact, is modeled by elements with reduced effective modulli. The stress distribution near 
the flawed region shows a tendency toward stress concentration, whose severity is related 



to the values selected for the effective modulli, which in turn can be related to the degree of burning 
applied. 

I would also like to inform you of some interesting and somewhat related work involved with 
the characterization and nondestructive evaluation of metal matrix composites. My colleague, 
Professor Blucher, has developed a new technique to manufacture metal matrix composites by contin- 
uous casting. We are interested in understanding the effects of manufacturing parameters, such as 
temperature, infiltration temperature, cooling rate, and matrix materials on the mechanical properties 
of these composites. We are also interested in predicting the mechanical properties of our composite 
using acousto ultrasonic technique and acoustic emission. We have acquired an imaging system and 
have used ICEPACK software to develop pattern recognition parameters to distinguish good and bad 
composites. We have broken some of our composites and have collected their acoustic emission 
activities. The acoustic emission activities of an alumina fiber reinforced 606 1 aluminum is shown in 
Figure 1 1. The results show that the activities of the specimens increase exponentially up to .22% 
strain. However, beyond .22% strain the activities rate changes.and becomes a linear function of strain. 
This change of activities may be related to the changes in micro mechanisms of failure. We are 
currently investigating the damage mechanisms at different strain levels, in order to be able to explain 
the change in the energy release mechanisms. One may postulate that the early events are associated 
with matrix yielding and fiber matrix debonding, and the later stages are associated with fiber fracture 
and linkage of cavities. Several specimens will be taken to different strain levels and then unloaded 
prior to their failures. These specimens will then be sectioned polished and studied under both optical 
and scanning electron microscope in order to reveal damage mechanisms at different strain levels. 

We hope this gives some information about our progress. I will provide you with a 
comprehensive report on each project as soon as students finish their theses. In the mean time, may we 
take this opportunity to wish you a merry Christmas and a happy and healthy new year. 

Sincerely yours. 


Hamid N. Hashemi 


and 


John N. Rossettos 
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Fig. 1 Schematic diagram of the model for the peel stress analysis. 
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Fig. 2 Effect of symmetrical central void size on the peel stress distribution. 




Fig. 3 Finite element model for the wave propagation analysis. 
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Fig. 4 The received wave form at node 1060, due to a broad band input pulse with the center 
frequency of 2 MHz for a specimen with 52% of the over lap symmetrical central void. 
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Fig. 5 Effect of symmetrical central void on the peak amplitude of the received wave at the node 1060. 
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Fig. 6 Failure locus of a tubular bonded joint subjected to tension/ torsion. 
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Fig. 8 Variation of the tensile strength with the heat damage time of fiber glass epoxy composites. 
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Stress & Events vs. Strain in Nextel 440/606 1AI Composite 
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Fig. 1 1 Events versus applied strain in an alumina fiber reinforced aluminum metal matrix composite. 


